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UNTERSUCHUNGEN ÜBER DIE KONZENTRATIONSABHÄN- 
GIGKEIT DER OBERFLÄCHENSPANNUNG UND DER VISKO- 
SITÄT DER LÖSUNGEN VON STARK POLAREN SUBSTANZEN 

IN POLAREN UND UNPOLAREN LÖSUNGSMITTELN 


Von M. Pucnuarik 


Physikalisches Institut der Schlesischen Medizin-Akademie, Rokitnica 
(Eingegangen am 21. Oktober 1953) 


In der vorliegenden Arbeit wurde die Konzentrationsabhüngigkeit der Oberflächen- 
spannung und der Viskositát der Lósungen von stark polaren Substanzen untersucht. In 
fast allen Fällen kann die funktionelle Abhängigkeit der Oberflächenspannung von der 
Konzentration durch die Szyszkowskische Gleichung dargestellt werden. Nur im Falle der 
Aethylalkohol-Benzol und der Pyrocatechin-Wasser-Lösungen ist die betreffende 
Abhängigkeit von ganz anderer Art.— Bemerkenswert ist auch das Verhalten der Parachore 
der Lösungen, die sich als streng additive Funktionen der Parachore der Komponenten 
nur für die Alkohol-Benzol Lösungen erweisen. < 


Zahlreiche Messungen an Lósungen von polaren Verbindungen in unpolaren 
Lósungsmitteln haben bekanntlich zum Schluss geführt, dass die molare Polarisation 
gelöster polarer Substanzen eine funktionelle Abhängigkeit von der Konzentration 
aufweist. Nur in verhältnismässig seltenen Fällen, wenn die gelöste Substanz keiner 
Assoziation unterliegt, ist die Molpolarisation konzentrationsunabhängig. In allen 
anderen Füllen nimmt entweder die molare Polarisation mit wachsender Konzen- 
tration ab, oder sie zeigt bei einem bestimmten Werte der Konzentration ein mehr 
oder weniger ausgeprägtes Maximum (Debye 1928). 
| Für eine tiefere Kenntnis der Natur der Lósungen ist es wichtig den eventuellen 
Einfluss der Assoziation auf andere physikalische Konstanten, besonders auf die 
Oberflächenspannung und die Viskosität, eingehend zu untersuchen und die Resultate 
mit denen der Polaritätsmessungen zu vergleichen. Das bisher vorliegende Unter- 
suchungsmaterial betreffend die Oberflächenspannung von Lösungen ist sicherlich 
unzureichend, es trägt einen zufallsartigen Charakter und braucht also einer Ver- 
vollständigung, Verfeinerung und Systematisierung, zumal die Messungen ver- 
schiedener Forscher oft nicht unerheblich voneinander abweichen. Im Zusammen- 
hang mit dieser Arbeit schien es dem Verfasser wichtig die Brauchbarkeit T klas- 
sischen Steighóhenmethode zur Bestimmung der Oberflüchenspannung von Lósungen, 
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deren Anwendbarkeit zu diesem Zwecke seitens mancher Forscher bestritten wird 
(Rehbinder 1932), näher zu prüfen. In nahem Zusammehange damit steht auch das 


Problem der Parachoradditivität von Lósurgen. 


Messungsmethode 


Die Bestimmung der Oberflächenspannung wurde nach der Steighöhenmethode 
der Flüssigkeit in einer Kapillare ausgeführt. Als Kapillare wurde eine bakteriolo- 
gische Pipette verwendet, deren Radius sich als ausreichend gleichmässig erwies. Die 
Steighóhe wurde mit einem Komparator auf 0,1 mm genau abgelesen, was für die 
Oberflächenspannungsbestimmung einen kleineren Fehler als + 0,6% ergibt. Die 
Viskosität wurde mit einem Höpplerschen Viskosimeter für Industriezwecke mit 
einer Genauikeit von 0,5% gemessen. 

Als Lösungsmittel dienten Benzol und destilliertes Wasser. Das „Benzolum purum 
frei von Thiophen“ der Industrie wurde nach einem Verfahren gereinigt und ge- 
trocknet, das bereits in einer früheren Arbeit des Verfassers verwendet wurde (Pu- 
chalik 1932). Die physikalischen Konstanten des Präparats waren folgende: Dielektri- 
zitätskonstante: 2, 282 + 0,001; Brechungsindex: 1,501 + 0,0002; Erstarrungspunkt: 
(5,5 + 0,1)°C; Oberflächenspannung: (28,9 + 0,1) Dyn/cm. 

Das Wasser wurde vor jeder Messung zweimal in einem Destillationsapparat aus 
Jenaer Glas destilliert und gleich nachher verwendet. 

Alle Messungen wurden ausgeführt in einem Ultrathermostat nach Höppler bei 
der Standardttemperatur von (20 + 0,01)°C. Bei Messungen an leicht oxydierbaren 
Lösungen trachtete man danach die Messungen möglichst schnell zu vollbringen, 
ehe eine merkliche Oxydation eintreten konnte. 


Untersuchungsmaterial 


Es wurden folgende Lósungen untersucht und die Resultate mit den Ergebnissen 
der in den Klammern angegebenen Forscher verglichen: 

Aethylaether in Benzol (Hertzen 1902, Hammick 1929), 

Nitrobenzol in Benzol (Hammick 1929), 

Aethylalkohol in Benzol (Hammick 1929, Starobinietz 1951), 

Aethylalkohol in Wasser (Traube 1881, Firth 1920, Shefield 1925, Buttler 1932, 
Guggenheim 1933, Gugel 1941, Zhukhovickij 1944, Cherneva 1947, Teitelbaum 1951), 

Phenol in Wasser (Goard 1925, Harkins 1925, Biegelow 1928, Swearingen 1928, 
Wagner 1929), 

Pyrocatechin in Wasser (Harkins 1925, Swearingen 1928), 

Resorcin in Wasser (Harkins 1925, Biegelow 1928, Swearingen 1928), 

Hydrochinon in Wasser (Harkins 1925, Swearingen 1928). 

Die Substanzen, die nicht in geniigendem Grade der Reinheit zur Verfiigung 
standen wurden sorgfältig gereinigt. Das Aether sulphuricus pro narcosi wurde z. B. 
nach Trocknung mit Kalziumchlorid mit metallischen Natrium behandelt und dann 
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destilliert. Die Oberflächenspannung des so 
+ 0,05) Dyn/cm. 


Der Aethylalkohol wurde vermittels vom Kalziumcarbid entwässert und dann 
destilliert, die Bestimmung des 


gewonnenen Präparats betrug (17,0 + 


kritischen Punktes des Gemisches von geprüften 
Alkohol und Petroleum gab für das Wassergehalt einen kleineren Wert als 0,1%. 


Die Resultate der Messungen sind in den Tabellen 1 bis 8 und in Abb. 1 zusammen- 
gestelt. Es wurden folgende Bezeichnungen angewandt: 
f; Molbruch der gelósten Substanz, 
fı Molbruch des Lösungsmittels, 
c Konzentration in Gewichtsprozenten, — ^ 
c,, Konzentration in Mol /Liter, ze 
M, Molargewicht der gelösten Substanz, 3 
M, Molargewicht des Lösungsmittels, 34 
с Oberflichenspannung der Lösung 32 
in Dyn/cm. 30 
d Dichte der Lösung in g/cm?, 28 
P Parachor der Lösung 
P, Parachor der gelösten Substanz, 
P, Parachor des Lösungsmittels, 
n Viskosität der Lösung in Centipoisen 
Die Resultate der Messungen für die 
Lösungen von Aethylaether in Benzol, 
Nitrobenzol in Benzol Aethylalkohol in 7 
Benzol wurden an der XIII Tagung der Mo 01 02 02 04 05 06 07 08 Q9 105 
Polnischen Physiker in. Kraków im De- i Fig. ] 
zember 1950 vorgestellt (Puchalik 1951). TE 
Gleichzeitig beschüftigten sich mit dem Problem der Konzentuationsabhängigkeit 
der Oberflächenspannung der Benzollösungen von primären Alkoholen die sovie- 
tischen Forscher Starobinietz und Rizhikova (1951). 


(C,Hj,0 — С.Н, += 20°C + 0,01°C 


Р, oe Pofa 


206,8 
207,2 
207,9 
299,1 
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Tab. 2 


C,H,NO, — C,H, += 20°С + 0 01°С 


Ри; + Pofa 

0,0000 

0,1821 30,0 216,2 

0,2545 . 31,1 220,3 

0,3935 32,8 228,2 

0,6074 36,0 240,2 

0,7603 39,3 247,8 

1,0000 42,9 

Tab. 3 


C,H,OH — C,H, += 20°C + 0,01°C 


Pıfı + Ра 


0,0000 


0,1537 193,8 
0,2389 187,1 
0,4506 170,4 
0,6136 157,4 
0,7496 146,7 
0,8811 136,3 
1,0000 
Tab. 4 


0,7894 
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Tab. 5 


C,H,OH — H,O t= 20°C + 0,01°C 


m — С,Н,(ОН), — H,O t = 20°C + 0,01°C 


Р, + Ри» 


Der Vergleich der Resultate des Verfassers mit denen anderer Forscher welche 


ң ISCH b 
ä sbestimmung benutzten, zeigt im All 
ee ess rim mit: Die Abweichungen sind jedenfalls 


sen der einzelnen oben erwähnten Forscher 


gemeinen eine genügende Ubereinst 
nicht grösser als zwischen den Ergebnis 
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untereinander. Damit scheint die Anwendung der Steighöhenmethode zur Ober- 
flächenspannungsbestimmung von Lösungen berechtigt zu sein. 

In fast allen untersuchten Fällen sind die о = f (c) Kurven, welche Ше funktio- 
nelle Abhängigkeit der Oberflächenspannung von der Konzentration darstellen, von 
normaler Gestalt, d. h. sie kónnen durch die bekannte Szyszkowskische Gleichung! 
o = og — b In (1 + ac) erfasst werden. Nur im Falle von Aethylalkohol-Benzol-und- 
Pyrocatechin-Wasser-Lósungen sehen die o = f(c) Kurven ganz anders aus. Die 
Kurve für die erstgennante Lösung lässt sich gut durch eine gerade Linie approxi- 
mieren und bei der gegenwärtigen Genauigkeit der Messungen kann man nicht 
bestimmt entscheiden, ob die vorhandenen Abweichungen von der Linearitát wirklich 
bestehen oder der mangelhaften Reproduzierbarkeit der Oberflächenspannugswerte 
zuzuschreiben sind. 

L. K. Staveley und Betty Spice (1951, 1952) ziehen auf Grund eigener Unter- 
suchungen über die Volumkontraktion von Aethylalkolhol-Benzol Lósungen, sowie 
auf Grund der Hoeckerschen Messungen der Dielektrizitätskonstante (Hoecker 1936) 
und der Bestimmungen der Gefrierpunktserniedrigung (Beckmann 1888, Peterson 
und Rodebusch 1926) einen wichtigen Schluss über die Assoziation der Alkohol- 
moleküle in der genannten Lösung. Danach sind die Alkoholmoleküle in ganz ver- 
dünnten Lösungen nicht merklich assoziert, in mehr konzentrierten Lösungen da- 
gegen bilden sie durch Wasserstoffbindung polymolekulare Komplexe in der Art 
von Dimeren, Trimeren und Tetrameren. Es liegt der Gedanke nahe in diesem Um- 
stand die Ursache der Abweichung der Kurve für Aethylalkohol — Benzol Lösungen 
von der gewöhnlichen Form zu suchen. Diese Deutung haben auch die sovietischen 
Forscher (Starobinietz et al. 1951) im grunde angenommen. | 

Bemerkenswerterweise zeigt in diesem Falle die у = f(c) Kurve ein allmähliches 
Ansteigen der Viskositätswerte bei zunehmender Alkoholkonzentration (Golik und 
Rabinowich 1950). Im Falle von Alkohol-Wasser Gemischen besitzt dagegen die 
n = f(c) Kurve ein ausgeprägtes Maximum, während die Oberflächenspannungs- 
werte allmählich mit zunehmender Konzentration abnehmen. In dieser Lösung tritt 
bekanntlich die Assoziation der Alkoholmoleküle gegen . die Wechselwirkung von 
Alkohol und Wassermoleküle zurück. 

Die endgültige Entscheidung, ob man die gewonnenen Resultate zur Unterschei- 
dung der verschiedenen Assoziationsarten verwerten kann bleibt künftigen Untersu- 
chungen vorbehalten. 

Die oben gegebene Deutung für die Alkohol- Benzol Lösungen versagt im 
Falle von Pyrokatechin - Wasser Lösungen. Kryoskopische Bestimmungen zeigen, 
dass keine nennenswerte Assoziation in dem untersuchten Bereich der Konzentra- 
tionen stattfindet. Eigentümlich ist auch das Verhalten der Parachore der unter- 
suchten Lösungen. Nur im Falle der Aethylalkohol- Benzol Lösungen erweisen 
sich die Parachore als streng additive Funktionen der Komponenten. 


20 — Oberflächenspannung des Lósungsmittels, со — Oberflächenspannung der Lösung, а, b = 
Konstanten. | 
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KPATKOE СОДЕРЖАНИЕ 


М. Пухалик, Исследование зависимости поверхностно1о натяжения и вяз- 


кости от концентрации растворов сильно полярных веществ в полярных и не- 
полярных растворителях. 


Была исследована зависимость поверхностного натяжения и вязкости рас- 
творов сильно полярных веществ от концентрации. Почти во всех случаях функ- 
циональная зависимость поверхностного натяжения может быть выражена урав- 
нением Шишковского. Только в случае водных растворов этилового спирта, бен- 
зола, пирокатехина соответствующая зависимость иного рода. 


Парахоры растворов только для алькоголь-бензолевых растровов оказались 
сильно аддитивными функциями парахоров составляющих. 
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SCATTERING OF 80° K—NEUTRONS 
BY CH;0H AND СН; MOLECULES 


By J. A. JANIK 


Physical Institute, Jagellonian University, Kraków 


(received January 19, 1954) 


Experiments on the scattering of slow neutrons by rigid molecules of a liquid show 
that the theory of Sachs and Teller may be applied to many liquids. The cross-section 
for the scattering of slow neutrons is then independent of the energy of the neutrons. 
This was verified for the СН, — molecule using 300°K — and 80°K - neutrons. Measure- 
ments of the cross-section for СН;ОН- molecules give the value (180 + 6) - 107?*cm? 
for 300°K — neutrons and (227 + 10) · 10 24cm? for 80°K — neutrons. The cross-section 
obtained from the theory of Sachs and Teller is 247 - 10-24 cm?. These results are inter- 
preted as the influence of torsional oscillations between OH and CH; groups in the CH,OH- 
molecule in agreement with the theory given by W. Kolos. 


1. Introduction 


It is well-known that the cross-section for scattering of slow neutrons by mole- 
cules depends on the geometrical structure of these molecules (Sachs and Teller 1941). 
The cross-section for scattering of a slow neutron by an atom bound in a molecule 
is according to the theory of Sachs and Teller 


д = og (илиьизи)* (То + Tr +.) (1) 


where o, is the cross-section for scattering of a slow neutron by the given atom in the 
case of zero-energy neutrons when the atom is rigidly bound to an infinitely great 
mass (in the case of protons g,, = 40, where о. = 20. 10-24 cm? is the cross-section 
of a free proton), То, Ті, ... are terms which depend on the energy of the neutrons 


: e n 1 
and on the energy of the translational motion of the molecules, и = 3 (uy + Me + из) 


and иу, из, Из are figures constructed in the following manner: 
A so-called „mass-tensor“ М of a molecule is introduced. It is defined by: 


557 r T 1 | (MY Se 
(М к = BT, , ij T 
ordinates of the considered atom as referred to the principal 


(167) 


ў 
where гу, rg; гз are the co 
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axes of inertia of the molecule, J}, Ip, Ig are the moments of inertia of the whole 
molecule as referred to these axes, М) is the mass of the molecule. 

One introduces а dimensionless tensor n = тМ-1, where m is the mass of the 
proton. The characteristic values of this tensor are n4, п», ng. Now 


E A Е ав 
DA pes BE ^ EE EAS E DS 


The theory of Sachs and Teller is valid for molecular gases and in the case when 
the neutron energy is low as compared with the vibrational quantum of atoms in the 
molecule but great as compared with the rotational quantum of the whole molecule 
(Eros < E, « [ЖЭР 

In a previous work (Janik 1953) it was suggested by the author that formula (1) 
may be applied also to many liquids. In this case it is necessary to replace the bracket 
expression (Т + T, +...) by 1. This corresponds to neglecting the energy of the 
translational motion of the molecules in comparison with the energy of the neutrons. 
The agreement of the experimentally obtained cross-sections with the theoretically 
evaluated from the theory of Sachs and Teller was demonstrated on the examples of 
several liquids (H,, CH; I, H,O). | 

According to our present knowledge of the structure of liquids (Frenkel 1946) this 
corresponds probably to the following picture: the above liquids are (at room temper- 
ature) in a solidlike state as regards translation, but in a gas-like state as regards 
rotation. They have then an almost free rotation possibility and a strongly hindered 
possibility of translation. 

The simplified form of formula (1) valid for liquids is 


a = O,, (и Ug ug и). (2) 


From this formula it follows that the cross-section for the scattering of slow neutrons 
by the molecules of a liquid is independent of the energy E, of the neutrons within 
some limits E, < E, < Е, (where E, < E, < E, < E, jj). 

The theory of Sachs and Teller as applied to gases or to liquids is valid only 
rigid molecules and perhaps also to molecules with free internal rotation (Janik 1952). 
It cannot be, therefore, applied in general to molecules with internal hindered rotation 
of one atom group with respect to another. 

Such a molecule is for exemple the CH,OH molecule. It is well known from 
spectroscopical data that the proton of the OH group in this molecule can rotate 
around the CO axis. This rotation, however, is not free, but it is rather a torsional 
oscillation on account of the relatively high potential barriers formed by the three 
protons of the CH, group. The frequency of this torsional oscillation (known from 
spectroscopical measurements) is about 250 cm7! (Borden and Barker 1938, Koehler 
and Dennison 1940, Burkhard and Dennison 1951). The energy of this oscillation 
is about 0.03 e V, i. e. of the same order of magnitude as the energy of thermal neutr- 
ons at a temperature of 300? K. 
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On account of this resonance it seems interesting to investigate the influence 
of the torsional oscillation on the scattering of 300° K-neutrons. According to the 
analogy between torsional and valence oscillations one may expect a lowering of the 
cross-section of the СН,ОН molecule in comparison with the cross- 
from formula (2) (Bethe 1937). 

Experiments on the scattering of 300° K-neutrons by CHOH molecules were 
made by the author (Janik 1953). The following results were obtained: the theoreti- 
cally evaluated cross-section of the CHOH molecule (considered as a rigid one) from 
formula (2) is 247. 10-24 cm?. The experimentally obtained value is (180 + 7). 

10-4 cm?. For the similar but completely rigid molecule CH,I the theoretically 
evaluated cross-section is 195. 10-24 cm? and the experimentally obtained one is 
(195 + 10) · 10-24 cm?. The lowering of the cross-section for methyl alcohol was 
interpreted as the influence of hindered rotation in the CH4OH molecule. 

The correct theoretical explanation of lowering of the cross-section was given by 
Kolos (1953). He obtained the following formula for the cross-section of an atom 
bound in a molecule with the internal hindered rotation 


section evaluated 


Cry mti = PO” Он, (3) 
т denotes here the set of rotational quantum numbers т, т’, j and о is the quantum 
number of hindered rotation. The indices 7 and f correspond respectively to the ini- 
tial and final states of the system. 6 рт; is {һе cross-section of an atom bound to 
a rigid molecule and 52 is a coefficient < 1. 

Kolos evaluated that for collisions with neutrons of 300°K (these collisions are 
elastic with respect to the hindered rotation, 
connected with great changes 
of t, for instance 0— 10,2— 11 
etc.), 52 is about 0,7 for a 
proton of the CH, group and 
near 0 for the proton of the 
OH group in the CH,OH mo- 
lecule. The theory of Kolos is 
then able to explain the lowe- 
ring of the cross-section. This 
explanation is only semiquan- 
titative as formula (3) is not Fig. 1 
equivalent to formula (2). А is 
very difficult averaging over the Maxwellian energy distribution of the neutrons and 
over the distributions of the energy states of molecules of the liquid is necessary. 

It follows from this theory that for slower neutrons (for example of 80°К) this 
lowering should be smaller, because the slower neutrons can cause only transitions 
connected with a small change of т. Then, the coefficient b? would not differ conside- 


rably from 1. 
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Taking into account all the above considerations it seems to be interesting to 
measure the cross-section of the CHOH molecule for 80°K-neutrons, as well as the 
cross-section of the CH,I molecules for 80°K-neutrons to verify the suggestion of 


its independence of the energy of the neutrons. 


2. Apparatus 


A source of photoneutrons from the Be (y, n) 2He reaction was used. 200 mg of 
radium in platinum tubes and needles were enclosed in a little cylindrical box filled 


. dh 
dE 


1 
N 


*400V 


tothe 
lonisation 
chamber 
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Fig. 2 Fig. 3 


with 0,5 kg of beryllium powder. Neutrons were slowed down to thermal energies in 
a layer of paraffin wax 6 cm thick (Fig. 1). The 80°K-neutrons were obteined in an 
arrangement similar to that described by Westcott and Niewodniczaüski (1935). 


In one side of the paraffin wax cylinder an opening was made containing a soft glass 
Dewar vessel. This vessel was filled with paraffin wax cooled by liquid air. From this 
cooled paraffin wax thermal neutrons of Maxwellian velocity distribution correspond- 
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ing to the temparature of 80° K were obtained. It was not possible to use common 
Dewar vessels, because they are made of hard glass containing boron and are not 
transparent for thermal neutrons. The energy distributions for neutrons of 80°K and 
300°K are given in Fig. 2. 

The neutron detector was a baron lined ionization chamber already described 
in a previous paper (Janik 1951) (Fig. 1). The ionization chamber was put within 
a cadmium box, which possessed in the front part a cadmium tube directed towards 
the Dewar vessel. In such a way only neutrons from the Dewar vessel could be counted 
by the ionization chamber. The external electrode of the ionization chamber was on 
the potential of 400: V. On account of the high per cent of dampness in the 
laboratory room. it was necessary to heat the inside of the cadmium box by an 25 
W bulb. 

The pulses from the ionization chamber were not amplified by an electrometric 
valve but they were conducted directly te the grid of the 6J7 valve of the first 
stage of the amplifier. The outline of the first stage of the amplifier is given in Fig. 
3. Two other stages are quite similar to the first one. 

The counting device contained an univibrator cooperating with the counting 
valve KT66. 

Fig. 4 contains sketches of the neutron source and the cadmium box containing 
the ionization chamber. On the first picture the neutron source is shown without 
the upper paraffin wax layer and the cadmium box without the lid for showing the 
beryllium powder and the ionization chamber. 

The scatterers (chemically pure CHOH and СН) were placed in flat glass 
vessels. The thicknesses of the liquid layers were: 0.156 cm, 0.232 cm and 0.261 cm. 


3. Results 


Cross-sections were evaluated from the measurements of the decrease of neutron 
flux intensity after the passage of a given layer of the scatterer. For each of two in- 
vestigated substances three values of the cross-section were obtained, together with 
the corresponding average errors. They correspondend to three thicknesses of the 


TABLE 1 


Cross-section per 


molecule for 300? Cross-sections per 


Cross-section per Cross-section per 


Substance | K-neutrons experi- | molecule for 300° | molecule for 80? | molecule theoretically 
mentally obteined in| K-neutrons experi- K-neutrons experi- | evaluated from for- 
the previous work | mentally obtained in | mentally obtained in mula (2) 


(Janik 1953) this work this work 


(195 + 10) 107?*cm? | (204 + 10) 107?*cm? | (193 + 12) 10-24ста? 


CH,OH 247 - 10-24 cm? 
3 


195 · 10-24ст? 
Сну 
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scattering layers. From these three values the weighted averages together with the 
average and probable errors for each substance were evaluated. 

The measurements for 300°K-neutrons were repeated and the new measure- 
ments for 80°K-neutrons were made. 

The results obtained are given in table 1. The errors given are the probable 
errors of the measured cross-sections. 


4. Discussion 


The results for CH,OH and СН, for 300? K-neutrons are in agreement with those 
obtained in the previous work (Janik 1953). 

The cross-section value for CH,I molecules for 80? K-neutrons is the same within 
the error limits as for 300? K-neutrons and in agreement with the theoretically evalu- 
ated one from formula (2). This result confirms the hipothesis of the independence 
of the cross-section of the molecules of a liquid of the energy of the neutrons. 

The results for CH,OH show a strong dependence of the cross-section of CH,OH 
molecules on the energy of the neutrons. Comparison with the theory of Kolos gives 
the following results: the mean energy of the 80° K-neutrons is about 4 times lower 
than the mean energy of the 300° K-neutrons. Consequently the 80° K-neutrons 
produce transitions connected with a small change of the quantum number т (0 — 4, 
1- 4 etc). For these transitions the values of the coefficient b are about 0.95 for 
the protons of the СН» group and about 0.4 for the proton of the OH group. The 
lowering of the cross-section in comparison with the rigid molecule will then be of 
a few per cent only. It can be seen from table 1 that the experimentally obtained 
result for the 80° K-neutrons agrees well with the theory of Kolos. 

The results obtained confirm the following details of the structure of the CH,OH 
and СН molecules: 


1. The geometrical models of the СН] and CHOH molecules in agreement 
with the spectroscopical data, 


2. The torsional oscillation of the OH group as referred to the CH, group in 
СН$ОН molecule. The frequency of this oscillation is about 250 cm-! and it is also 
in agreement with the spectroscopical data. 

The distributions of the energy of the neutrons at both temperatures are very 
diffuse and therefore it is not possible to evaluate from these experiments the height 
of the barrier of the hindering potential in the СНОН molecule. 

It seems interesting to make similar experiments using monoenergetic neutrons 
obtained in a neutron velocity selector. 


My thanks are due to the Ministerstwo Zdrowia PRL as well as to Dr E. Wyrobek 
for permitting me to use 200 mg Ra belonging to the Instytut Onkologiczny in Kraków. 
I wish also to express my gratitude to Professor H. Niewodniczaüski for his kind interest 
in this work. 
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KPATKOE СОДЕРЖАНИЕ 


M. Яник, 


Pacceanue нейтронов температуры 80° K na молекулах СНОН 
и CH. 


В одной из предыдущих работ была поставлена гипотеза, что теория Сакса 
и Теллера эффекта рассеяния медленных нейтронов на твёрдых молекулах газа 
относится тоже к многим жидкостям. В этой работе гипотеза эта была прове- 
рена для молекул CHsI употребляя нейтроны температуры 80° К и 300° К. Теория 
Сакса и Теллера не применяется однако к молекуле СНОН по причине внутрен- 
ной заторможенной ротации групп OH и СН; по отношению к себе. Полученные 
результаты на эффективное сечение для CH3OH согласны с результатами опу- 
бликованными в предыдущей работе и с поданной Колосом теорией влияния 
заторможенной ротации на рассеяние медленных нейтронов. 
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QUENCHING OF PHOTOLUMINESCENCE OF SOLUTIONS 


Ву A. JABŁOŃSKI 
Physics Department, Nicholas Copernicus University, Torun 
(received January 25, 1954) 


The present theory of quenching of photoluminescence of solutions is based on 
following assumptions. The quenching molecules in solutions carry out constantly irre- 
gular oscillations about their actual equilibrium position, which they change from time 
to time. The frequency of changes of the equilibrium position depends ceteris paribus 
on the viscosity of the solution and of its temperature. To every distance of the equili- 
brium position of a quencher from that of the luminescent molecule belongs a certain 
time proportional quenching probability. The above distance of the equilibrium posi- 
tions is assumed to be discrete — the quencher may be present in the first, second, and 
so on, shell constituted of the molecules of the solvent surrounding the excited lumines- 
cent molecule (the „shell model“ of the luminescent centre). The behaviour of a system 
of such luminescent centres is described by a system of differential equations. The theory 
is applied to some simple cases, and expressions are obtained describing the decay of the 
total (i. e. emitted in all directions) photoluminescence intensity as well as expressions 
giving the quantum yield of photoluminescence as a function of concentration of quen- 
chers and other factors. Apart from quenching by quenchers. also the „inner quenching“ 
of luminescent molecules and the quenching during the „initial shock“ (disturbance in 
the motion of nuclei in the centre caused by electronic transition and, lasting a very short 
time) are taken into account in these expressions. The problem in question is solved 
rather generally for the case of solid solutions in which the equilibrium positions of the 
quenchers do not change. Only approximate expressions are obtained for the more 
involved case of liquid solutions. Several expressions given so far by different writers 
result as particular approximations from the present theory. 


1. Introduction 


The quenching of fluorescence of gases is described adequately by the well 
known Stern-Volmer (1919) formula: ? 
1, 0 

N иттен ө с бр | 9 

(c. f. also P. Pringsheim 1949, p. 90), where I, and v, are respectively the intensity of 

fluorescence and the mean life time of the excited molecules in absence of quenching 

and z, the number of quenching collisions per second, 2,70 being proportional to the 

_ pressure p of the gas causing the quenching. The mean life time of the excited states 

| .(175) 
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1 
of quenched molecules is given by— = — + 2,, whence 
T To 


ee (2) 


Eq. (1) and (2) show that 
т. @) 


The success of Stern-Volmer’s theory induced several writers (Vavilov 1928, 1936; 
Vavilov and Frank 1931; Sveshnikoff 1935 and the present writer in his doctor 
thesis (unpublished)) to look for an explanation of queching phenomena in liquids 
in the same way, the only difference being the replacing of the relative gas kinetic 
velocity of the molecules in gases by the relative velocity of Brownian motion of 
fluorescent molecules and quenching molecules („quenchers“) in liquids. A formalism 
developed by Smoluchowski (1917) in his theory of coagulation was videly used in 
these theories to describe the quenching collision rate (theories of diffusion con- 
trolled quenching). 

The deviations of experimental results from those predicted by the above theories 
induced Vavilov to introduce a „configurational quenching“ in addition to the „col- 
lision quenching". His expression for the relative yield 

e 9c 


0, си (4) 


where c is the concentration of quenchers апа о and k certain constants, seems to 
agree well with his experimental results. w has the meaning of the effective volume 
of the excited molecule and 


12 _ RT (o + с») (c1 + 95) 


37 01 05 


, (5) 


where 9, and 9, are the „kinetic“ radii of the excited quenching molecule and o the 
effective radius of the excited molecule. According to Pringsheim (1949, p. 326) 
the assumption from which (5) is derived is valid only in a restricted sense. 

If only the static (configuration) quenching is taken into account there results 
from (4) the equation | À 

0 = Qe, (6) 
which is identical with F. Perrin’s (1924) expression for efficiency of photolumines- 
cence versus own concentration of the solution (case of self quenching). A more com- 
plicated expression for this last case was given by Vavilov. 

According to Rabinowitsch (1937) collisions in liquids occur in sets, called co- 
ordinations or encounters (the so called cage effect). A molecule has a reduced pro- 
bability of collision with a second molecule after a collision with a first one, but has 
an enhanced probability of hitting the same molecule several times (cf. Williamson 
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and La Mer 1948 and Pringsheim 1949, p. 326). Clearly, the situation in liquids differs 
fundamentally from that in gases. In gases the configuration of atoms or molecules 
changes so rapidly that it is possible to consider collisions as independent events 
(no Wahrscheinlichkeitsnachwirkung!), whereas it is not so in the case of liquids. 

The purpose of this paper is to develop a theory of quenching and of related 
phenomena in which the above property of liquids is taken into account and which 
could embrace all different cases of quenching including even that of quenching in 
gases (as a limiting case). This theory may seem to be somewhat formalistic since it 
is developed without any special assumptions concerning the specific mechanism of 
quenching. It seems, however, that such a formal theory, free of any uncertain assump- 
tions, may prove useful. It shold be noted that polimerization, which plays in certain 
cases a paramount role in self-quenching of photoluminescence (Levshin 1927, 1934), 
is not accounted for in this theory, and, if the theory is to be applied to such cases, 


the influence of polimerization (say, dimerization) on the yield of photoluminescence 
has to be taken separately into account. 


8 2. Shell model of luminescent centres 


The theory, which has to be developed below, is based essentially on the fol- 
lowing model of the luminescent centre in solutions. We assume that a luminescent 
centre in a solution consists of an excited luminescent molecule enveloped by shells 
consisting of monomolecular layers of molecules of the solvent. Among them some 
quenchers may be present!. Different centres in the same photoluminescent solu- 
tion differ, in general, by the occupancy of iheir shells by quenchers. Thus, there 
are always different kinds of centres in a solution in which quenchers are present. 
If the solution is a solid one, the quenchers remain infinitely long in their original 
shells. In liquid solutions, however, there exists a finite time-proportional probability 
that a quencher may pass from one shell to another. This probability depends on the 
viscosity of the solution and on its temperature. We assume, besides, a time-propor- 
tional probability of quenching the excited molecule by a quencher, depending, 
ceteris paribus, on the shell in which the quencher is situated. А quencher situated 
in the shell nearest to the luminescent molecule possesses, of course, the largest proba- 
bility of quenching, a quencher in the second nearest shell has a smaller one (but 
not necessarily nil) and so on. We thus assume that quenching in solutions Eon due 
to collisions of quenchers with excited molecules, but rather to „encounters during 
which the quencher remains a certain time in a particular shell .of a centre. During 
this time the quencher oscillates around its equilibrium position (cage effect). The 
above assumption of the time proportional probability of qnexehing for a centre with 
a quencher in a particular shell (which is, strictly speeking, the mean value of that 


1 Only quenching phenomena will be considered in this paper. The writer intends to treat depolar- 
ization of photoluminescence by means of the same model in a subsequent paper. 
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N N = 


probability) can be valid only if the time of stay of the quencher in a shell is large 
compared to the period of oscillations. 

We now consider a luminescent solution as a system constituted of N excited 
luminescent centres immersed in a solvent, which acts also as a heat bath. We divide 
all the centres in groups according to the occupancy of their shells by quenchers: 
a group of centers without quenchers, a group with one single quencher in the first 
shell and so on. If there occurs a transition of a quencher from one shell to another 
in a particular centre, this centre changes his group. We may call this change of 
group „transition“ of a centre from one group to another. In the state of statistical 
equilibrium there are as many direct as reverse transitions between any two groups 
per second. If the distribution of centres between groups differs at a certain time from _ 
that of equilibrium, the number of transitions in one direction becomes unequal to 
that in the reverse one. This unequality lasts until the state of equilibrium is reached 
by the system. 

Let N; be the number of centres belonging to group i and n the number of groups 
in which all the N centres are divided. Then 


Y N, — N. | (7) 


We assume the time dependence of N; to be given by 
Ме 2 (4; Nj — AjN) — (y + Fo + V;) N; + Ay i=123..n, (8) 


where 4; ed 4; denote respectively the „transition probability“ of a transition of a 
centre d SUE j to group ? and that of the reverse transition, y the 
transition probability of electronic transition giving rise to the photoluminescence 
band, W, the „inner“ quenching probability (i. e. the probability of quenching in- 
dependent of the presence of quenchers in the centre), W, the probability of quen- 
ching by quenchers characteristic for group i, and A, the number of new centres 
belonging to group produced by excitation during one second. If the last two terms 
in (8) would be equal to zero, A, №; — 4,; N; = 0 in the state of equilibrium of the 
system. Thus 


en (9) 
Apo NU 


where № and № correspond to the equilibrium distribution. Hence? = Ay 


* Eq. (8) appears to be somewhat similar to that used by Cox (1950) for description of irrever- 
sible changes in an ensemble of systems (cf. also Thomsen 1953). There exists, howerer, a substantial 
difference between them. While eq. (8) has to describe the behaviour of a system, that of Cox describes 
the behaviour of an ensemble of systems. The present formalism seems to fit better to our special pur- 

` pose, and enables us to introduce easily terms by means of which the emission and quenching of photo- 
luminescence аз well as excitation of centres can be taken into account. It may be noted that, while 
Àj = A,jin the formalism of Cox, Aj; + Àj; in general in the present theory. 
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The general theory will be illustrated below by some simple applications. Let 
us consider a system excited instantly at the moment t — 0 (the case in which the 
duration of the excitation is very short compared with the time of decay of the excited 
centres) in which no further centres are produced at t >> 0. In order to describe this 
case we must put A, = 0 in (8). Let N, be the number of (excited) centres at t — 0, 


and let us denote N;/N, by p,. Thus » р; =latt = 0, but M р;< lattz 0, because of 


the decay of centres due to the emission and the quenching of photoluminescence. 


For A; = 0 and N, = p;N, the system of eq. (8) becomes: 
р;= (Aj Pj — дур) ОМ + ТР) р» і = 1, 2,... п. (10) 


$ 3. Simple applications of the general theory 


The simplest possible case (the crudest approximation) is that, in which centres 
are considered as luminescent molecules with one shell only?, in which there may 
be at most one quencher present. In this case all the centres may be divided in two 
groups only: the centres with one single quencher in its shell forming, say, the group 1, 
and those with no quenchers the group 2. The set of eq. (10) reduces to 


pi (y + Fo +F + Àj) P1 — Арз = 0, 
Рә T (y + Wo + А) Pa — 4р = 0, (11) 


where A, denotes 4,5 and A, denotes 121. W is the quenching probability by one 
quencher present in the shell. 
The solution of (10) is of the form 


-@+ Бе, 


TY e 1 
ра = a4 € € + by е0. (12) 
where 
Я-А W 
amy Р A uL 
and A 
AcTAY,F(u—A) (EF (13) 
B GS ANE 
Substituting (12) into (11) and taking (13) into account we find 
ау ры "i À1— А5 d. vs | (14а) 
ENT | 359 ai | 


3 Instead of the notion of a „shell“ the notion of an „effective volume“ should ا‎ d 
this approximation, since it may often prove to be much larger that the actual volume of а monomole- 


cular shell. ; 
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b 1 14-42 W 
eee] did 
If for t =0 
Ра E SG 
Pi Po 
then аз + ba = C (a4 + b1) = 
Fu F ply 2: E ira ا‎ 
a? [e 2 ФЕ 
Непсе 


желу; Ay—A И 
[| И | (15) 
By means of (14) and (15) the first eq. (12) becomes 


er дүе Ай ор DELS 
м (5 - aes Е в, | e-(a+Bl pee] e d 
(16) 


E 17] 
e a ha +. в) sal ы 


Hence the second eq. (12) becomes 


е) арен 
" 18 
desee] ® 


The constants C and A have to be fitted to the initial conditions. Let v cm? be the 
volume of the shell, п cm“? the concentration: of the quenchers and 


» — nv «&l. (19) 
If the system is in a state of statistical equilibrium at t = 0, the probability* of finding 


a centre without quenchers at that time is psy = e7”, and that of finding a centre 


with one single quencher руу = ve^" (we neglect for the sake of simplicity the centres 
with more than one quencher). Thus 


ee (20) 


ты Pa at bee RZ 


4 Cf. Smoluchowski 1904 or 1907. 
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Since pp = ve^" = А 2f, 


ДЇ єз 2g ` (21) 
From (9) and (20) we have A, /A, = 1/ v. We introduce 


À can be considered as a constant only for a given system at a given temperature. . 
Generally it depends, ceteris paribus, on concentration of the quenchers and on 


viscosity and temperature of the solution (cf. next section). 
Eq. (17) and (18) become finally 


E I; EIE ? =| eet + L e = еті, (23) 


| 


2p 
1 1 
| = (1-1) +» (1+2) 
Pa = > МЮ YN e +8 + |1 сы ج‎ е-@—#‹| (24) 


where a nad f are given by (12), which by means of (22) become 


À 1 
e-1 2 (1+1) + 


ANCE oe 


From (23) and (24) we obtain an expression describing the decay of the system, and 
thus also the decay of the total photoluminescence,5 i. e. the decay of the total radiation 
emitted by the system in all directions per second, if the excitation of the system 
occured instantly at % = 0. The expression reads . 


| А? (1 — ») + (24 +4) 
e | en are ae И ШОП 
( | ч (25) 
(1—») Y--212--»40-— | 
+ (1+7 + 7j М —! 


p is the probability of finding a centre in the excited state at time 2 after the moment 
at which excitation has taken place. 
We shall now consider some particular cases of (25). 
If W = 0 (no quenchers in the solution) eq. (25) leads to 
р =e” (1+ 0) е9 + wi д (1 — »2) e-& +, (26) 
оиа o a a ————————— 
5 The decay of particular components of polarized photoluminescence will be treated by use of the 


same model in a subsequent paper. 
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The factor е-” (1 + v) ~ 1 — »? instead of 1 arises from neglecting groups of centres 
with more than one quencher in the shell. The decay in this case is purely exponential. 


If A — 0 (the case of a solid solution) 


D = e” (уе-@ + We W)t е0 +5). (27) ` 


This equation shows that if there are quenchers in a solid (or very viscous) solution, 
the decay of photoluminescence does not proceed according to a simple exponential 
law, even in the simplest cases in which there are no depolarizing factors, and in 
which the assumption of centres with one single shell is sufficient for an adequate 
description of the behaviour of the solution. In this particular case (A = 0 — solid 
solution) one can easily solve the system of eq. (10) in which centres with more than 
single quenchers are also taken into account. The system for this case becomes 


р; =— (y + Fo + Wi), і =0,1,2... (10а) 
where i denotes the number of quenchers in a centre. The probability- p;, of finding 


i quenchers in a centre at t = 0 is e” ә! (cf. Smoluchewski 1904 and 1907.) 
The solution of (10a) reads 


1 Е 
р; = ре tet MIE = е" = соте оа (27a) 
VW; denotes here the probability of quenching of an excited centre with i quenchers 
(possibly И; == iW, at least in some cases). For i =0 W,=0. The decay 
curve (27a) is even more complicated as that given by (27). It approaches an expo- 
nential curve for t — оо. 
For W— œ eq. (25) becomes 


PIER THES v. 12 as (28) 


Independently of the value of A the decay proceeds in this case according to the 
exponential law. If A = 0 (solid solution) the quenching is purely configurational: 
all centres which possess quenchers are quenched immediately after the excitation 
has taken place. 


If (28) is applied to the photoluminescence of gases % => *) has the mea- 
v 


ning of the number of quenching collisions per second. 


$4. Quantum yield of photoluminescence 


Apart from time proportional quenching probabilities, which were assumed in the 
foregoing sections, there exist certainly quenching processes, the probability of which 
can by no means be considered as time proportional. Every electronic transition in 
a luminescent centre leads the centre into a state of nuclear motion in general very 


different from that of thermal equilibrium. (cf. Jablonski 1931) The thermal ваа. 
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brium is reestablished very quickly, but before it is reached, there is more (sometimes, 


much more) opportunity of quenching than in the state of reestablished equilibrium. 
We should like to call this kind of quenching „quenching by initial chock“. Let us 
denote by Q, the quantum yield of photoluminescence, which would be observed 
if no quenching would taken place except that due to the initial chock. The quan- 
tum yield Q with all quenching agents taken into account can be calculated as follows 


© = Qo y | pdt, (29) 
0 


where y denotes the probability of the electronic transition responsible for emission 


of the photoluminescence band and p = 3 Р; (сї. eq. (10)) is the probability of finding 


а centre in the excited state at a time t after the excitation has taken place (р = 
= l for = 0). For our simplified, probably even oversimplified, one shell model, р 
is given by eq. (25). It has the form 

р = Ae (FPN |. noch (30) 
where A and В are constants the meaning of which may be found by comparison 
with eq. (25). a and f are defined by eq. (13a). From (29) and (30), we get 


which is the general expression for the quantum yield for the one shell model. 


Instead of discussing this general expression, we shall rather consider some simpler 
particular cases. \ 
For the case of a solid solution (A = 0) eq. (27а) and (29) lead to 


N. к ада 3 
It may be remembered that » = vn (cf. (19)), where v is the volume of the shell and n 


the concentration of the quenchers. eu. 
If only the first two terms are taken into account, eg. (32) can be rewritten as 


e ln FEW (83) 
ЖОҒА АЦ y Н уы сы 

; 1+» УР 
where О, denotes the yield for v=o. Eq. (33) is equivalent to expression (4) 
given by Vavilov. 


If Wz y+ Wo, (33) becomes approximately 

Q ~= Qoe”, (34) 
Perrin’s expression for the case of self-quenching. 
is the concentration of luminescent molecules 


which is equivalent to eq. (6), 
y in this case denotes Nv, where N ; К 
and v as before the volume of the shell (ог rather that of the „effective volume“). 
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Since even eq. (32) must be considered as approximate (because it is derived 
for the one shell model of the centre), eq. (33) and (34) are a fortiori only approximate. 

Let us now consider the case in which A = 0, but W— оо and restrict ourselves 
to a two group approximation only, centres without any quenchers forming one 
group and centres with опе or more quenchers. forming the other one. Only centres 
without quenchers can emit photoluminescence; p for this case is given by eq. (28). 
From (28) and (29) we obtain 


у 


E77 AT (35) 
Queis Fr. ( a8 3 


Q = Ое" 


À 1 
As was mentionned above, one must терага іп this case 2 1+ Se the number 


of quenching collisions per second (A is certainly a function of concentration of the 
quenchers and of the temperature of the solution). Remembering that у = vn and 


1 | 
putting A + 1) = kn, ме obtain from (35) 
y 


алау m 
О = Qe y + Fo + kn ` (36) 


For vn <1, (36) becomes approximately 


en (87) 
where the notations are self-explanatory. Eq. (37) is equivalent to eq. (1), i. e. to the 
well-known Stern-Volmer formula (from which also eq. (2) and (3) follow). 


85. Concluding remarks 


The main feature of the present theory is, apart from the introduction of the 
»Shell model“ of the luminescent centre, the assumption of a time proportional quen- 
ching probability depending on the shell in which the quencher is situated. The case 
of solid solutions were considered more exactly than that of liquid solutions. Expres- 
sion (27a) derived for the one-shell model for solid solutions can be easily generalized 
by making use of the many-shell model®. The case of liquid solutions (A = 0) is much 


$ Let po (пу, п»... п„) be the probability of finding a centre with n, quenchers in the first shell, 
п» in the second and so on, and n, in the last shell z at t = 0, and let v; be the volume of the shell i. 
If the concentration of the quenchers is п cm "3, the distribution at t = О is: 


n! > 5 
Ро (п1, пз... п.) = GE EEE ЕТЕГІ TEE ng 1-- Up По 


5 че (38) 
па n4! (5-2) ! k=1 i=1 ‚ | | 


j=1 
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more involved. It is described by the system of eq. (10). In order to make clear how 
the formalism works in this case it seemed expedient to ayoid tedious calculations. 
Accordingly, we restricted ourselves to the consideration of the simplest possible 
case, that of centres with one single shell only (which is, of course, equivalent to the 
familiar „effective volume“). Moreover, only centres with no quenchers and centres 
with at most one quencher were taken into account. The system of eq. (10) was thus 
reduced to the system (11) of two simultaneous equations. Their solutions were dis- 
cussed. One may hope that even so crude solutions may prove sufficient for an appro- 
ximate description of the phenomena in question, but there is no difficulty of reaching 
better approximations also in this case (i. e. the case of 2+ 0). It may be of interest 
to investigate the nature of A from the point of view of kinetic theory. Such investi- 
gation is not made in this paper. 

It follows from considerations given in $3 that the decay of photoluminescence 
does not follow in general a simple exponential law, even if the total intensity (not 
that of a single componenet?) is measured. The experimental values of the mean 
duration of photoluminescence 7 may depend on experimental method used? and 
may differ considerably from that defined by 


oo 
аура 
0 


The method of treating the quenching phenomena given in the present paper may 
seem somewhat formalistic. It seems, however, to fit well to the properties.of liquids. 
[t does not need any assumptions on the specific mechanism of quenching. 


If we assume the probability of quenching of a centre with shells occupied by quenchers as given by 


z . . . 
(38) to be УР, where ІР, is the probablitity of quenching by n; quenchers situated in shell i, the decay 
1 t 


ісі 
function becomes 


x А 
Фа Уһ (my, пз... Mz) exp | — (+ Wo+ у) Wi) ; (9) 
D 


Nis п... ї= 1 
and the yield 
У 
Q= 93 Ро (пл, п... п) = o ЕИ (40) 
Nyy Ne.. Nz yt Wot EV, 
i=l 


Presumably there will be no need of using such involved expressions as (39) and (40), and eq. (27a) 
will prove sufficient for all practical purposes. ге 
7 The decay of single components. of polarized photoluminescence was already treated (Jablonski 
1935, 1936). It will be treated again by use of the shell model in the next paper. 
8 The methods in which whole decay curves are obtained, as that used e. g. by Franck and Prins- 
heim (1943), are clearly superior to those in which only the mean duration of роон 18 


measured. j 


- 
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Severäl expressions, given by different writers, for the yield of photoluminescence 
as a function of concetration of the quenchers, are obtained here as particular appro- 


ximations from the present theory. 


KPATKOE СОДЕРЖАНИЕ 


А. Яблонский, Тушение фотолюминесценции растворов. 


Некоторые явления, связанные с тушением фотолюминесценции растворов, 
можно описать формально системой дифференциональных уравнений. Рассматри- 
ваются некоторые простые настные случаи. Из решений системы для этих слу- 
наев получены приближенные выражения для затухания и для выхода фото- 
люминесценции, из которых некоторые, более упрощенные, совпадают с известны- 
ми выражениями предложенными другими авторами. 
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ON THE MOTION OF BODIES 
IN GENERAL RELATIVITY THEORY 


Ву Г. InreLo 
Institute of Theoretical physics, University of Warsaw, Warsaw 


(received February 19, 1954) 


The object of this paper is to give a simple derivation of the non - Newtonian equa- 
tions of motion following from the field equations in general relativity theory. This is 
achieved through the use of Dirac's Ó functions in the energy-momentum tensor. All 
calculations (with the exception of the most trivial ones) are included. 


1. Introduction 


The problem of motion in general relativity theory seems to be sufficiently 
clarified. The method used originally by Einstein, Infeld and Hoffmann (1938) 
consisted in treating matter as singularities of the field. Thus the field ‘equations 

Rag = 0 
prescribed the motion of singularities, provided the new approximation method was 
used. Almost simultaneously Fock (1939) dealt with the same problem of deducing 
the motion of gravitational bodies from the field equations and in his paper obtained 
Newtonian motion; later the non-Newtonian approximation was also obtained by his 
method (Petrova 1949) and agreed with the result of Einstein, Infeld and Hoffmann 
(1938) and Einstein and Infeld (1949). a 6 

The comparison of Fock’s and our method shows one essential similarity. 
In both methods the same approximation procedure is used. But.the two methods 
also show essential differences. The first consists in the fact that whereas we use the 
field equations for empty space, regarding masses as singularities, Fock uses the 
gravitational equations with the energy-momentum tensor. This difference I regard 
as unimportant ideologicaly for the following reason. It is jt that no charge or 
gravitational body can be properly described by a singularity. We can use the singu- 


larity picture only if we describe the field outside the particle, where the energy- 


momentum tensor vanishes, It is obvious from our derivation that it would have been 


almost unchanged had we used surfaces covering the regions in which the energy 
momentum tensor does not vanish. Indeed Fock's method gives precisely the same 
result as ours though the mathematical technique is very different. 

$ (187) 
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The other difference between Fock’s method and ours is more essential. It con- 
cerns the choice of the coordinate system. To compare Fock’s paper and ours on 
this score we must distinguish between two versions of our paper. In the first version 
(Einstein, Infeld and Hoffmann 1938) a coordinate-system was used similar to that 
employed by Fock. In the second version (Einstein and Infeld 1949) apparently no 
coordinate system was used at all, yet the equations of motion were exactly those 
obtained before. Yet Fock regards the choice of a coordinate system as essential. The 
difficulty was clarified later (Infeld and Scheidegger 1951, Infeld 1953). Indeed the 
approximation method itself forms implicitly a stringent enough coordinate condition 
to determine the equations of motion up to the sixth order. 

Our derivation of the equations of motion required the application of many 
ideas, the use of the „new approximation method“ being only one of them. The 
technical part turned out to be quite difficult so that the calculations were not given 
fully either in the first version or in the second one.I do not believe that the derivation 
of Fock and his collaborators is simpler than ours. After a long familiarity with the 
problem, after having had the pleasure of discussing it with Professor Fock, I was 
convinced that there must be a short cut to the derivation of the equations of motion. 
Indeed the present paper is — certainly from the technical point of view — a great 
simplification. The proof of this lies in the fact that all the ideas and all the calculations 
(leaving out only trivial steps) are contained in one, not very long, paper. In some 
respects this simplification was achieved by a synthesis of the ideas of Fock — 
Papapetrou (1951) and our own. 

I shall explain this by an example. In the classical theory of gravitation we usually 
distinguish between two kind of equations; either 


А ф = 4ло | (1,1) 
if matter is represented by a continous distribution, or 


if matter is represented by singularities. In (1,2) we are usually looking for a solution 
with central symmetry and singular, say, for г = 0. In this case the solution is: 


Da 
p e r (1,3) 


I could say, roughly speaking, that Fock’s method corresponds to the assumption 
of (1,1) whereas our papers assume (1,2). Since we do not know the real distribution, 
neither of these methods properly depicts reality. The use of (1,2) means: if the two 
bodies are a great distance apart so that central symmetry of the field can be assumed 
approximately near one body, then the exact knowledge of the density distribution 
is not essential and equation (1,2) can be used. | 

Nowadays more and more often one writes (1,2) when thinking about its solution 
(1,3) in the form 
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where ó(s() is the 3-dimensional Dirac ô function, zero outside the point of origin 
and such that the volume integral over the whole space is one: 


ЕСЕ dax = 1 


The use of Dirac’s д — function is not of great advantage here. Yet one sees 
that (1,1) and (1,2a) cease to be two quite different cases. By spreading out the ó 
function one obtains a continous distribution. Taking о in (1,1) greater and greater 
and keeping the total mass constant, one arrives at the 6 function in the limit. Mathema- 
ticians object to the way physicists use these functions but since Schwarz's work it 
is known that a rigorous theory of these functions exists. 

Thus the use of б — functions forms a bridge between Fock's method and ours. 
But from the example given it would seem that the use of ó functions can give no great 
simplification. True enough, we can generally write as a solution of (1,1) 


(М) Ex — == dix 
5-2 


апа applying this to (1,2a), we have: 


9 x' m N 7 r 
ф = —т о 28 $4 = = | да (х) Чад = =——. 
|к— «| [x | | | í 


But this can hardly be regarded as an essential simplification. Yet a skilful use of ô 
functions in General Relativity Theory simplifies immensely the derivation of the 
equations of motion. 

Let us return for a moment to the classical equation (1,2) or (1,2а) and assume 
that we have s particles, the masses of which are: 


1 2 5 
UOS 1s, eso M 


and their motions are described by 


1267 s 2) Were 
Ek Е... Е where k =1, 2, З and &* = £ (t). The solution of (1,2) is: 


p=1 


о Po — B. 


Now to obtain the same solution from (1,2a) we have to introduce on the right 


hand side of it 
t, vp 
4л > т бз) 


p=1 
where 


E р р TM. 
бз. = д (x1 — 1) à ) س‎ £?) ò (x — 63). 
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This notation will be used throughout the paper. For the sake of simplicity we 
shall assume s = 2, that is the two body problem throughout. In this case the Newto- 


nian equations of motion are: 


17 152 a y: 

-— 2 E (1,34) 
r 

22 21 2 Y. 

m PER е) (1,35) 
Me 


1 і 2 

г? = (EF = 2%) (Gi zu £^ 
Yet our aim is to deduce these and then more exact equations of motion from 
field equations alone. Therefore it may be useful to remark that we can write even 
the Newtonian equations of motion — with the use б functions — is a „field form“, 


Indeed instead of (1,3) we can write: 


291 17772 
131 15921203 TOME Ver Ek в рю 8 
оь О ee =O) а^ rae 


r3 73 


Integrating this equation over a volume {2 surrounding the first singularity we obtain 

(1,3a) and integrating over a volume (2 surrounding the second singularity we obtain 

(1,3b). Finally integrating (1,4) over a volume containing both singularities we obtain 
11 2 2 


mtm =0 


that is, the conservation laws for linear momenta. 


2. Notations: the gravitational equations. 


Since in the greater part of our work we shall have to separate space and time, 
our notation will not be the usual four — dimensional one. We introduce the con- 
vention: Latin indices take the values 1, 2, 3, and they refer to space coordinates only. 
Greek indices refer to both space and time, running over the values 0, 1, 2, 3. Repe- 
tition of indices implies summation. 

The expression 


баға etc. stands for сх еїс. (2,1) 


At infinity the gravitational field takes the Galilean values N, that is: 


| Uist а Din 0; Nog = 1 (2,2) 
We write: 
Em = 1 bh g^ ot №" (2,3) 
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where ,, represents the deviation of space — time from flat space and is not assumed 
to be small. 


The A" can be calculated as functions of h,, by means of the relation 
Em 8^ = д, (2,4 


The gravitational equations (with the use of the energy momentum tensor) are 
in natural units: 


1 
Ras сұл сые ӛлі, (2,5) 
ог: 
1 
Rag = —8л er 9 бав г) (2,6) 
We shall use tensor densities: 
Ru; ar V-g Ra, ; Тов т. V-g Tap (2,7) 
and the equations that we use here will be: 
1 
Rog = — 8л (r.— 3 Sap т) (2,8) 


Now the left hand side of (2,8) is well known. We шау recall: 


л 
- e {co = 


The question is what shall we use for the right hand of (2,8)—that is, for the 
energy momentum tensor? It is here that Dirac’s ö functions comes in. 

We shall assume that the right hand side is always zero with the exception of 
singular world lines of the s partieles. That is, we shall assume: 


Sp 
T = у Т" 
= | (2,10) 


PP Р P РРР р. 


T" = m(t) $; Tom = mi"; T" = т" д 


That is, in a coordinate system in which the p — th particle is at rest, the only com- 
р 


р 
ponent of the energy momentum tensor density is T® = тд. We assume, that the 
world-lines of the particles never intersect. These assumptions are essentially the 
same that we have made previously (1938, 1949) although they were not explicitly | 
stated, because the energy momentum tensor was not used. But the solutions consi- 
dered previously in the lowest approximation were those of the field with singular lines 


(2,10). j 
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3. Equations of motion. 


Writing: 
1 
Rt — E g” В = G^? (3,1) 
we know that in consequence of the Bianchi identities we always have 
from which follows: 
Ta. Bd 0 (3,3) 


In our case, in which 1“ is given by (2,10) we can write (3,3): 


Р 
Y r5,-0 (3,4) 


Let us now take the integral of (3,4) over the region sourrounding the p-th singula- 


rity. We have then: 2 
[т dix = 0 (3,5) 


ы 
9 


These are 4s equations. They are ordinary differential equations, that must be satisfied, 
if the field equations are satisfied; They are a consequence of the field equations. 
The only argument appearing in these equations is x° = t. Both, T"? and the gravi- 


РР 
tational field depend on m, &* and their time derivatives. Thus we have extracted from 


the field equations (which are partial differential equations) 4s ordinary differential 
p p 
equations which must be satisfied by the functions: m, & p = 1,... s. 


4. The approximation method. 


The problem before us is to solve our field equations and to deduce the equations 
of motion. This we shall do by a new approximation procedure. Let us assume a func- 
tion p (x", A) developed into a power series in the parameter A (for small A): 

9 (4,2) = 2° p+ Mod opt... = Е Мо (4,1) 
X 0 1 2 1 

The indices below indicate the order (/ in A! is always the exponent, not theindex). 

If the function ф varies quickly in space, but slowly with x^ then we are justified 
in not treating all its derivatives in the same fashion. The derivatives with respect to x? 
will be of a higher order than space derivatives. We can formalise the procedure by 
introducing an auxiliary time т, 

TEA (4,2) 
so that derivatives with respect to т can be treated on the same footing as the space 
derivatives: 

дф дф | 


Pike о э А = ^9» (4,3) 
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We conclude: the „stroke differentiation“ of a quantity with respect to x can 
be replaced by the „coma differentiation“ with respect to т if the power of A with 


which this quantity is associated is simultaneously raised by one. To express this 
explicitly we use numbers under zeros, written after the comma, 


21 س‎ 21-1 
5 h mnjo = hs Or: Zen =, 42427 (4,4) 
21 ДӘРІ 21 21 


тп,00 
2 


From now on, all differentiations will he with respect to (т, x", x, x3) and they will 
be denoted by commas: 


ШЕК ie mte TR (4,5). 


1 


Thus we shall develop all functions that appear in the field equations in power 
series in 4. We start with the Ts’. Because of (2,10) we have: 


Te = A2 T+ 22 TR. 
2 4 

To" = Аз Тот .. 45 Тот... (4,6) 
3 5 


у prio 


ла Т" ҒАУ +... 
4 6 
By convention we assume 


p p Р 
т = 22 т + Ат |... (4,7) 
2 4 5 


Р p ; 
that is, we start T% with T° Then T^?" = 7992” Therefore Т" starts with 7°”, and 
2 3 


P р 2 ы . . 
T"" with Т since T"" = T° 2" E". In the power series development we take into 
4 . 
account only odd or only even powers od A (The reason for doing this was explained 


in the paper of Infeld 1938). 
Even to start with m of the order two is not purely a convention. Indeed in the 
2 
Newtonian approximation we have: 


mass X mass 


mass X acceleration = А 
ГА 


Since the acceleration is of order two, the orders of both sides will be equal only if 


the order of mass is also two. 
The linear terms in R,, аге: 


1 ; 
Roo ze E hoojss ES hosjos Sr 2 hss,00 E L 00 


T2 1 Ji l ; 
Ron = T€ Roniss SP 2 Роз E 2 "зов "g 2 hano + L 9n 
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1 1 1 1 
Вт = — 2 hmnjss + 2 Amsins ЕТ: 9 Рта — 9 Рап 
1 1 1 
+ x mie — у! — 5 Promo 
1 : 
zi m Роот tr L mn (4,8) 


where the L’,,in (4.,8) denote the non-linear terms. 

Because of (2,6) and (3,6) the lowest order of Ку and Кии is also two whereas Rom 
begins with three. The same applies to hg, h,,,, and hg, which begins with the order 
two, and three respectively. 

All the functions that will appear later are obtained from the h’s by summation, 
multiplication, differentiation. Thus to every component, the following rule applies 


throughout: Апу component having an om number of zero suffixes will have 


odd КЕ қ 
only powers of A in its expansion. 
even 


9. The Newtonian equations of motion. 


We shall not now proceed with the general theory, but instead we shall try to 
find the equations of motion in the lowest (Newtonian) approximation. Later the 
non-Newtonian approximation will follow. These two examples and the final remarks 
should make the general theory quite clear. We start with the equation 


1 1 РР 
жы p Lue mm) en бе. E 
R 9 100 8л (r^ 2Т ) 5 epi 
bp (5,1) 
ом = 87 m д 
The solution of this equation is: 
t b 
hy = — Уу 2mri (5.2) 
2 p=1 
P / 
where r is defined by (2,10). We shall also write briefly: 


and in the two body case, which we shall assume for the sake of simplicity: 


ф=Л-+ 8 (5,4) 
i2 22 
f=—2mr-!; g=—2mr} 
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It is also worth remarking, that according to (2,3) and (2,4), we have 


100 — — эс йе» 
; 200 ф (5,5) 


But let us mention also that so far there is no reason to assume that m is constant. 
2 
Indeed this will follow from the equations of motion. These equations are, for the 


first particle: 
1 ; 1 
f pesto pev = е Mr а = 
I, Fir |) E (5,6) 
о 


where i = О апа T" = 0 are: 
0 п $ | m $ 
m +H TS. БД Т"" + 2 2 T" + lol О (5,7a) 


Р: тт +21 redi rm = 0 (5,7) 


Now (5,7а) is at least of the order four. However (5,7b) is at least of order three. 
The integral of Т” т gives nothing since it can be changed into a surface integral 
over a surface on which 7” vanishes. Thus what remains is: 


1-1 1 
f (m ) dax = то = 0 (5,8) 
2 2 % 2 р 


1 2 
We see therefore that m and т аге constant. The theory does not tell us whether m 
2 2 2 
should be positive or negative. Because there is attraction, and not repulsion between 


p 
gravitational masses, we take m to be positive. Indeed by taking the first particle and 


2 
1 


removing all others, we see that m is its gravitational mass, since for large r the 


2 
1 


field is that of a particle with gravitational mass m. This is the same constant of inte- 
2 


gration that appears in the Schwarzschild solution, since our field for one particle 


is that of a Schwarzschild singularity when r is large. 
Let us now calculate the three equations of motion from (5, Та) and 5,6) up to 


the fourth order. We consider the two body problem and denote: 
2 
Ө = nee: (5,9) 
The equation (5,7a) written out explicitly up to the fourth order is: 


о То qe. 5,10 
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1 
The integral over the region <2 of the first expression is zero. The integral of the second 


expression is: 


1 1 
(mn = > (5,11) 
Now as to the last integral: 
1 1 1 1 UT 
= | Ф,Т%4ух-- | fs ûm dax + 5 | 8, 0m dax (5,12) 
lI £ 2 Ја 2/71 
2 о 2 
The contribution of the first expression on the right hand side on (5,12) equals zero. 
Indeed the contributions to this integral from x* — 7* = a* and а^ — 7* = — ай 


are the same in absolute value but opposite in sign. This is an essential point, because 
Ж, becomes infinite for x* = n*. The integration of the second integral on the right 
hand side in (5,12) is performed by putting x = n* into g,. Denoting 


8,, = 8, for x = "n, 


2 
2m _ 9Е() 
о 2, = ( = )- pi (5,13) 
U) 
we have finally: 183 
long 1 т т т 
п 1а (27 614) 
ир 


that is, the Newtonian equations of motion. 
Let us summarise briefly what we have done up to now. From the field equa- 
tions we have found hop. The equations of motion of the third order gave us the in- 
2 
$ 


formation hat m is constant. The equations of motion of the fourth order gave us 
2 


the Newtonian equations of motion. Thus the character of our field singularities 
determines the equations of motion. 


6. Transition to the next approximation. 


We wish to find the next approximation. This means that we wish to find equa- 
tions of motion up to the sixth approximation. They are: 


1 1 
ip (ге, дий? т“) di x = 0 (6,1а) 
LOT 1 
| (т, 24. т") dis X = 0 (6,1b) 
< Q : 


The first equation will determine m, whereas the second will give us the equations 
4 
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of motion up to the sixth approximation, that is beyond the Newtonian approxima- 
tion. But to be able to do that, we must know the field that appears in (5,7), if taken 
up. to the sixth approximation. The field components that we must know are: 


hoo» lmn hom hoo (6,2) 
2-499] Б 
Р 1 
But before calculating these components ме сап evaluate m from (6,1a). This we shall 
4 
do now, turning to (5,7b). We have: 


0 1 1 2 2 
in Ta omo T^^ — o, (пат) тен i) (6,3) 
[Slr mare - Le (2321 
00,7, = | 0] 7 Аа A +m o) (6,4) 


We take the three dimensional volume integral of (6,3) over 2 The result is: 
1 


Be nu (6,5) 


Now we can make use of the Newtonian equations of motion, since the deviation 
from them will give a contribution to (6,5) of the seventh order. Thus we can write 


instead of (6,5) 
1 1 
—2nirir- — m (tip) 
2 2 


Now the integral of (6,4) can be evaluated in the following way: 


EN ires ET 
5 Ro di) X се) Noe XUI a EEO 
0 Q Q 


(6,6) 
De 1 VER: le Цеа M 
pac (em) + lem д.1 dax = y mE + т() 7) 
Q 
Therefore the integral over (5,7b) calculated up to the 5th order gives 
1 1 ады МЫ? ) 
-2 — _ т? + < = 0 6,7 
(2 бада тет аталы | (6,7) 
If we suppress the arbitrary constant that could be reabsorbed by m, we have: 
2 
12 
11 rd 1871 22 
1 ARE IRB (6,8) 
net 


Now we return to the problem of calculating (6,2). According to (4,8), (2,6) and 


(5,1) we have: 
1 1 1 1 ТӨ” 
— j a + rore іне — g hun р Pa — дере (69) 


ч 
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We are looking for a solution, which in the case of a one body problem would be 
spherically symmetric and Galilean at infinity. Such a solution is: 


(6,10) 
The next step is to calculate Ay. Again (4,8) and (2,8) gives: 
3 


1 1 1 1 ТЕЛ. 5442 
= hou + > Raum +5 һы — > has = Sr lmn"ö+merö) (6,11) 
3 3 42 1 rA EE 2 2 


- - 


Introducing (6,10) into (6,11) we have: 


1 1 1.0 che: ИБ 
— 5 hons + 5 По — Qo = —8л{тзү" 0+ тё" д (6,12) 
2 2 2 2 
The solution which we shall adopt here is: 


hon = —2fi" — 2g è" (6,13) 
3 


This is not the only possible solution. About a more general solution and its influence 
upon the equations of motion we shall say something in the last section. 
The calculation of Ло is more troublesome. We shall give it in the appendix. 


4 
Here we shall only quote the partial result, that is only those expressions appearing 
in Ақ that influence the equations of motion. These expressions аге: 


2 ка 
1 ні лин CELUM т 
he cy 8 — he За tz (6,14) 
r pd 


We shall now proceed to calculate: 


< 


7. The equations of motion up to the sixth order. 


These are: 


m n m m - 
f | T". + Т, + Hi +2 ia T” 4- К zl dix (7,1). 
о 
The first expression in (7,1) gives zero. Because of (6,8) the second, multiplied 
by 2-5, gives; 
i3 dua Арысы Gm Ms - 1 1 
-2 рт mm m — À- m m y iwi o к „ут 
m" mT aan ma der елаз var 272,077 (7,2) 
ihe даса 
ty mg, en 


Or, if we introduce the Newtonian equations of motion: 


: E 1 т 1 u.“ 1 1 ^ 4 
ж Too dos minui ni, mE Ein 0; 
33 


1:4 1 m 2 : 1 
= 2-2 ўт т ns — pm —. = ут 
зерен EEO), e 
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The third expression gives: 


| ee ix 
ТТ Шр ДЕС а-ы 
2 


nA = РАІ 
= D] i БЕЛДЕ "d 


The fourth expression gives: 


2 1 1 172 1 
—mn" я д (homn + Ian Фо опт) (зух = —4тт { — тр (ж 
ж 


о 
DT 1 2 (7,5 
+ (2) вё—2тт (1) oa" 4 

Г K T жт 


The fifth — and last expression—is slightly more complicated. First, we have there 
the old Newtonian expression: 


Т%-- 22 Т» (7,6) 
(n) 7+ | 


The first expression is well known; the second (integrated) cancels out with the 
second expression in (7,2), because of the Newtonian equations of motion. What 
remains now is to calculate 


Е а : h E BE 
i oof T” (3х = М 2 79 Зоо 2 5 (е 3% 
о 4 2 


1 1 Le ee : 
In the first expression 5 ? Pom only E Em Will give us significant contributions. 


(This fact is explained more fully in the appendix). homp gives only (—2g Eg p; and 
finally Ay gives only the expressions collected in (6,14). 
4 


Thus we find: 


Jem =4 C), emen C) 
1100); 2 2 \7/ yr" TI» (7,7) 
Q 4 
E24] 112 : 
esnin(1) „rrr rit 


Now we add (7,3), (7,4), (7,5), (7,7) and obtain non - Newtonian equations of mo- 


tion:! 


12 
1 We reabsorb the 775 by substituting new units for т and т, т: old т = A+ new т; old mass = 


= А-2, new mass. 


t 
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1 


2 
t үй ХІ sg Әже? LES m m 1 
fi o 5 e E a es] E 


T а Br mew 1 em (7,8) 
de non omes к Ne | T +5 Отит ê” 


The equations of motion for the other particle are obtained Ьу replacing 
12 225 
mmt by тті? 

respectively. 


8. The coordinate conditions. 


Let us introduce a transformation: 
0% = x? + 2۴-1 مه‎ 
Ecl (8,1) 
x* = oi + Aa, 
k 
It follows from this transformation: 


ж. 
hom = Pom + 40, 
EHE Е Ael 


hos = И ae Пыл 3 Qum (8,2) 
k k k 


ру аа 
hoo = Foo + 2 40,0 
k k k—1 1 


Thus the change of the coordinate system will influence the equations. 
Let us now introduce the changes: 


hom vr hom + т 
3 3 3 
leg > Ryan AE amn T anm (8,3) 

4 4 4 4 

hog hoo + 2 20,0 
4 4 31 


The only possibility of influencing the equations of motion through this new choice 
of the coordinate - system is of influencing the equation of the sixth order by the 
choice of a ар. Indeed a, enters into (7,5) but there the a,,„, and — a, „и cancel each 
other. The only other place in which a, enters is (7,7) where we have 


1 
— homo + 2 hoo, 


But here again the contributions from (8,3) cancel each other. Thus the addition of 
an arbitrary a, „to hon does not change the equations of motion. It introduces а change 
3 3 


in o and the two contributions cancel each other out. But the choice of a, and a 
4 т 
ва? са 4 
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will have an influence upon the equations up to the seventh and eighth order. There- 


fore it must be possible to choose the coordinate — system in such a way the equa- 
tions of motion given here will be the exact equations of motion. 


Appendix 


If we wish to calculate fo and avoid infinities, we must assume, instead of the 
4 
1 2 1 2 


Dirac functions б and д spread out densities ó(a), ó(a) depending on a parameter а 
1 2 


and going to û and ô if a— oo. Therefore instead of ф we shall assume g™ to be 
the solution of a differential equation: 


BR =\8 л (m 5(а) dq s Ха) (А 1) 
2 2 


In writing we shall suppress the a’ s but it must be understood that all quantities have 
been made finite by spreading out Dirac's д functions. The equations of motion do 
not have to depend on б.) that is on the internal structure of the bodies, as long as 
„spherical symmetry“ is assumed. 

Let us, therefore, turn toward the calculation of Rog. According to (2,9) we have 
up to the fourth order (with the absorbed As): ` 


1 1 3 1 1 
< hoos + Боо — g P00 + о ФФ — 3 ФФ (A 2) 


and because 


Ross =2 Фо» (A 3) 
3 


we have: 
1 1 1 1 h AM 
R® = (1 —39) Кю —— 3 Фа uy 949 FPP TZ Фо — ALS (А,4) 


Now the right hand side of the equation 


1 
RY = —8л (мт) (А 5) 
up to the fourth order is: 
1 1 lo. 
в [тетю ттт 


1 1 1 | 
DEN yo TO Ts 
са E 2 i 24 +3 4 
- Therefore the equation of the fourth order is, because of (6,7) and (2,10): | 


jud Ве ME 
hs = Фао 299a + Ps Pat 9n B T mo Le d 
4 y 


(Аб) 


(м апа М being the arbitrary constants missing in өт) 
4 4 
We remember that what we need from hoo» is the calculation of if Room m йс 
1 
Therefore the spherically symmetric expressions in hy do not xum any contribu- 
4 
tions. Thus we can ignore all the expressions with 5 on the right hand side of (A 6). 


Also the expressions f,f,or f, д. Therefore instead of (А 6) we can write: 


hoos ~ foo +80+ 2&5 & + 27.8 + Er&ır 


12 
2772 26373 
ЛЕЗДЕ 8+8 | 
4 
Let us start with the last expression 
252 
b„=8nMÖö 
4 4 
Then we have: 
2 
2M 
4 
а 
r 


and 


"e 
Lie) n 

where the proportionality factor is a constant. We know, however, that the T 
cient at Я” is the mass 7 of the body; therefore all expressions giving a proportiona- 
lity factor times 7” must cancel each other and the constant M was added to insure 


4 
this cancelation. Therefore we may ignore all expressions which give a constant x & i 


‚or a constant xj". We may ignore fio) because for а > oo 


om which follows 


Motion of Bodies in General Relativity Theory 203 


and 


1 54 1 1 ид 1 ys = nf 1 
2 Cm и (3% = — т 1 (To) m да зух = т ( 1 ) ЙЫ д асзух 
эт 


r 


which gives a + oo an expression proportional to 5* and has, 
Therefore: 


therefore, to be ignored. 


1 2 
: O12 2 
Рао» ~ Во + 2,6 + 2/ Bs + 8, Bis + 8л г тд — I i) 


2 
We omitted here the last expression 8л М 6 since its meaning is equivalent to 
4 


А А З 1 
that of ignoring allexpressions: constant х N =— Es constant X ё „. Finally the ex- 
pression 2/6, gives in the neighbourhood of the first singularity for а > oo: 
1 
—2 (7), go 


Therefore again the corresponding integral gives for a + оо constant x £, and 
can be ignored. Finally 


28,8 
requires a little more consideration. Let us write 


2 
Sr = 25, Si 8 m? ф(а) yha) ss 


where 
2 2 2 2 
ф(а) ss = — 4л8(а) = гі (гу) = ple!’ 
if 
2 
гр = & 
Therefore: 
rs 41| EE” dr — c| 
C 
Soc 8m? [Је] 
2 
Y C 
Im = ms FF a Г EE” dr + (5) | 
r эт. 
wheie: 


С = fe £^ dr 
0 


2 . LJ = 
so that at г = 0, s,,, would be finite, for a finite a. As an example let us take 
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Е-і--ет“4 
ЖАЗ; 
E 
r r 
2 
ЕЕ” Жеміс a3e-ra + а?е-?та 


Then: A ы 
fee ar =o; [гга --2--% 
r 0 


1 
In the limit for sufficiently great a we again have an expression (s = 
Г m 


which has to be ignored. 


— 22, if we add the expressions that can be ignored. 


Finally g,g, — 5 


Therefore: 


12,2 
1 325230277 та 
hoo,ss ~ 8:00 + (5#) +в termê 778) 
4 ЕС 


Therefore we have for h the expressions given in (6,14). 
4 


KPATKOE СОДЕРЖАНИЕ 


Л. Инфельд, О двиҗениц тел в общей теории относительности. 


Целью работы является вывод вне-нютоновских уравнений движения. Уда- 
лось, употребляя функцию д в тенсоре энергии-импульса, так упростить вывод 
этих уравнений с общей теории относительности, что все вычисления (кроме 
тривиальных упрощений) приведены в работе. 
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NOTE ON THE Z-DEPENDENCE 
OF THE CROSS-SECTION FOR PAIR PRODUCTION 


By А. Z. HRYNKIEWICZ 
Institute of Experimental Physics, Jagiellonian University, Kraków 
(received March 1, 1954) 


А new interpretation of experimental results of previous works of the author on this 
subject is given. The Z-dependence of the cross-section for pair production by у-тауз 
of ThD (2.62 MeV) was found to be consistent with the theory of Jaeger and Hulme. 


In previous papers (Hrynkiewicz 1952, 1953) experimentally obtained values 
were given of the ratios of the cross-sections for electron pair creation in Al, Cu, 
and Ag to the cross-section in Pb. These results were rather inconsistent with the 
values calculated on the basis of the theory of Jaeger and Hulme. This discrepancy 
was especially distinct іп the case of Al. А new interpretation of the results of the 
measurements is now given. i 

The measurements, as reported in the papers mentioned above, consisted in 
registering by two G.—M. counters the rates of coincidences N for aluminium, copper, 
silver and lead foils of various surface densities s, when irradiated by a beam of y — 
rays of ThD (2.62 MeV). The curves N versus s were drawn and the ratios of coin- 
cidences Np,/N as a function of s were evaluated for Al, Cu and Ag. These ratios 
were found to depend on the surface density owing to the different absorption of 
electrons in foils of the same surface density of different materials and on the scat- 
tering of Compton electrons. Undesirable coincidences due to the direct Compton 
effect were eliminated by a suitable geometry of the experiment. i 

As the number of the created electron pairs is directly proportional to the surface 
density of each foil and the effects due to the absorption and scattering of electrons 
are proportional to s?, the ratio of the rates of coincidences due to electron pairs 
could be obtained by extrapolation of the values of Np,/N to zero thickness of the 
foils. The thinnest foils for which measurements were performed had the surface 
density about 10 mg/cm?. Therefore, for the extrapolation of Ng/N the shape of 
each curve N(s) between 0 and 10 mg/cm? must be postulated. In the previous papers 
it was assumed that the ratio Np,/N = R(s) should approach a constant value as 
we pass to zero thiekness of the absorber. This means that the condition 


(205) 
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dR 
= o 


must be fulfilled. This condition suggested the manner in which the interpolation of 
the curves N(s) in the interval 0 — 10 mg/cm? was performed. However, as it is 


shown below, the condition (1) does not hold. 
ТҺе rates of coincidences for the two foils can be written in. the form 


'Then 


Coincidence rate counts/hour 


5 


0 10 20 30 40 
Surface density mg/cm? 


Fig. 1. Rates of coincidences as function of 
surface density of various absorbers. 


gives the ratio of the numbers of created electron pairs and the first derivative of 
В at s = 0 does not disappear 
dR ad — bc 
ЕНЕ E : we 


с 


That. means that the curve R(s) must reach the R axis at an angle different from 90°. 

Not being restricted by condition (1) the curves N(s) between 0 and 10 mg/cm? 
could be drawn as smooth extensions of their course for higher s. These curves are 
shown in Fig. 1. 
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Curves R(s) evaluated from curves N(s) are represented in Fi 
in Fig. 2. F 
values of R(0) = (N>,/N) P 2 1g rom the 


s= o the ratio of the cross-sections for pair production o/op, 
may be easily calculated. 


A 
zx 


Ratio of coincidences R(s)» 


——À 
0 10 20 30 40 
Surface density mg/cm* 


Fig. 2. Ratios Np;/N as functions of s for 
Al, Cu, Ag and Pb. 


In Table I the newly evaluated experimental results are compared with the ra- 
tios of the cross-sections calculated from the theories of Bethe & Heitler and Jaeger 
& Hulme. The limits of the errors of experimental values are now larger 
than those given in the previous papers since the interpolation of the curves N(s) 
between 0 and 10 mg/cm? without condition (1) is less definite. . 


Table 1 


ІІ” 


Ge cn 


Element : 
experimental 


Al 0.025 0.020 _ 0.021 + 0.004 

Cu 0.101 + 0.012 

ГЕ 0.267 + 0.025 
10 


——————-—_-_—_—_—————-—-—-—-—-С+ —— = St" = = = — 
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In the case of aluminium it is difficult to decide with which theory the experi- 
mental result is in better agreement. In the case of copper and silver, however, the 
results obtained confirm distincly the theory of Jaeger and Hulme. 

The results of the investigations of the Swiss group (Hahn et al. 1952) and Dayton 
(1953), published meanwhile, are also consistent with the theory of Jaeger and Hulme. 


KPATKOE СОДЕРЖАНИЕ 


A. Гринкевич, Заметка о зависимости от 2 зффективною сечения на обра- 
зование пар. 

Подана новая интерпретация экспериментальных результатов прежних работ 
автора на эту тему. Получена зависимость от Z эффективного сечения на обра- 
зование электронных пар ү-лучами ThD (2.62 MeV) согласная с теорией Erepa 
и Хольма. ; 
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ON THE MUTUAL INFLUENCE OF SPHERES IN VIBRATING AIR 


Ву W. STASZEWSKI 


Physical Laboratory, Maria Curie Skłodowska University, Lublin 
(received March 3, 1954) 


The force acting on small spheres in a Kundt tube placed at different distances 
from each other with their line of centres normal to the vibration vector has been investi- 
gated. When the spheres are close together this force is shown to be one of attraction 
or repulsion according to the size of the spheres as well as to the frequency and intensity 
of air vibrations. The contradiction of R. S. Cook’s and E. N. da C. Andrade’s results 
is thus explained, the conditions of their experiments differing considerably. 


Kundt’s experiments on dust figures in glass tubes have been known for nearly 
ninety years. Kundt (1866) found that with particles of dense material (sand, iron 
filings) sharp ridges were formed between the nodes in a glass tube when the air in the 
tube was excited to strong resonant vibrations. The explanation ОҒ: the dust 
tube phenomena was given by Kónig (1891) who calculated the forces acting on two 
spheres S, and S, situated in a vibrating non-viscous fluid. He found the following 
expressions for the components of the force acting on the sphere 5; 

SRO Me r. vs 


X, = и sin 9 (5 cos? 0 — 1), 


2r* 
ЕЕ 
2. 2 O ST Go: ~ 
: nS 
the Z axis being assumed parallel to the vibration vector (Fig.1). 
These formulae give ; 
Зло 113 га? vy? 
for д = 0° X, = 0; Z, = 21—20. 
Рз Волт fo n га vo” ; 2, = 0, 
2r* 
i. e. longitudinal repulsion and transversal attraction. For other values of д the 
direction of the force does not coincide with the line of centres of the spheres, becom- 
ing normal to it at 9, = 54240” and à, = 125^20' (Fig. 2). M 
The force is proportional to the volumes of the spheres and inversely proportion. 
to the fourth power of the distance between their centres. 
(209) 


and for 9 = 90? Жа = 
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König’s theory was considered satisfactory by several writers (е. g. Thomas 1913, 
Robinson 1913). Some authors however were not satisfied with his argument based 
on the assumption of vortex-free motion of the non-viscous fluid surrounding the 
spheres. The problem was taken up by Andrade in 
two important papers published in 1931 and 1932. 
Andrade excited sound waves using a loudspeaker 
fixed at, and closing one end of the glass tube, 
the loudspeaker being energized by a modern valve 
oscillator of acoustic frequency. He carefully inve- 
stigated the motion of solid particles in the vibrating 
air column at different frequencies and intensities. 
He discovered the phenomenon of the antinodal 
disc, described and explained many details of the 
formation and motion of ridges, but rejected Kónig's 

| theory as inadequate to explain the regular spacing 

Fig. 1: Two spheres in the sound field. ; H А : с 
ее тай of the ridges, the indeption of ridge formation when 
rı = radius of the sphere S the velocity of the vibrations exceeded a certain 
гә = radius of the sphere $, _ minimum, etc. Experimenting with very light spheres 

г = distance between the cen- about 2 mm in diameter prepared from pith Andrade 

узу UC э, observed that when a single row of spheres was cement- 

v — velocity amplitude of the А = 

ыан ИШ ed to the wall of the tube in a line normal to the 

9 = angle between the line of &Xis, the free spheres formed themselves into a row 
centres and the 7 axis. near the fixed row, and they returned to the same 

position after the tube has been tapped. This beha- 

viour of light spheres cannot be accounted for by the assumption of forces derived by 
König for vortex-free motion of the surrounding fluid. Andrade showed that suffi- 
ciently vigorous vorticesaround 
single particles and particle 
groups are essential for the for- 
mation of ridges. He made a 
qualitative study of the mutual 
influence. of two cylinders! 
suspended in the acoustic field 
and found that in the trans- 
verse position the cylinders 
attracted each other at small — a ae 


: © 
Q ©; 


105° 


o 
со 
< 


120° 


distances when the equilibrium р ig 2. The force acting on the sphere S, at different angles to 


position of the suspended cy- the line of centres with the vibration vector. 
linder was inside the vortex 


system initiated by the fixed cylinder in the absence of the second cylinder, and 


WEIT ڪڪ‎ ане ee ee о mm 


J) small spheres were also used. 
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repelled each other if the equilibrium position was just outside the vortex 
system. 

Earlier investigations of Cook (1902) who studied the force acting on sealing 
wax spheres, either lying on glass or suspended on spider’s thread, in the Kundt 
tube excited by a vibrating brass bar are in contradiction with Andrade's results. 
Cook found namely that when the spheres were situated transversely to the vibration 
vector there was repulsion at distances smaller than the radius of the spheres, and 
attraction at larger distances. Andrade suggested that this repulsion might be due to the 
electrification of the spheres, Cook's observation being that the effect was less notice- 
able (not more as stated by Andrade) with glass than with sealing wax spheres, but 
he did not pursue the matter further. 

The experiments described below were undertaken with the view of explaining 
the discrepancy between the results of Cook's and Andrade's investigations. 


APPARATUS 


To excite the tube into resonance a 6 watt electrodynamic loudspeaker L (Fig. 3) 
of Philips, type 9884, was used. It was energized by alternating current of any desired 
frequency from an oscillation generator O—G connected to the Philips 10 watt ampli- 
fier A,, type 2864. In most 
experiments the generator 
used was one made by T. 
Kopaczek, Warsaw (type 
P05). Some experiments 
were performed with a 
Philips oscillation genera- 
tor, type GM 2307. The loud- 
speaker was connected to 
the glass tube RE of 35—40 
mm in diameter by means 
of a large conical wooden 
tube C, a rubber sleeve | 
between the narrow end of Fig. 3. Experimental arrangement for studying the -mutual 
the cone and the glass tube influence of spheres in the acoustic field. 


providing a convenient air- T жі, 
tight connection. Ву the other end of the tube the Philips electrodynamic micro- 


phone M, type 9585, was placed. The microphone current was amplified by another 
valve amplifier A, and its wave form was examined on the screen of the cathode ray 
oscilloscope O. The approximate value of the frequency could be read off the dial on 
the face of the oscillator. The voltage on the loudspeaker was read Do the voltmeter Г. 

The objects to be observed in the acoustic field were placed insido the tube ЕЕ 
а few cm from the open,end E and lighted by a 30 watt lamp provided with a suitable 
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convex lens. The observations were made either visually or photographically, different 


optical arrangements being used. 


EXPERIMENTS 


I. Kundt dust figures in fine cork powder. 


At sufficiently high frequency and suitable intensity the ridges are stable and 
very sharp, being but one particle thick. As observed by Andrade »the individual 
particles can be plainly seen, well separated, in the ridges which are sharply in focus“, 
Two photographs of the ridges formed at frequency 560 c. p. s. (Figs. 4 and 5) confirm 
this observation. The separation of a large number of particles suggests strongly the 
presence of repulsive forces acting between the spheres situated transversely to the 


vibration vector- 


Fig. 4. Kundt dust figures in reflected light. Fig. 5. Kundt dust figures in transmitted light. 


II. Forces between spheres resting on glass in the Kundt tube. 


Small spheres made of Carnauba wax and of glass were used in these experiments. 
The wax spheres were prepared in the way described by Andrade. A conical vessel 5 cm 
high, opening upwards, was made of steel and provided with a hole about 0.25 mm 
in diameter. It was mounted on the axis of a cream separator filled partly with wax, 
heated electrically to about 100°C and, when the wax was melted, set in rapid rotation. ; 
The stream of molten wax split into perfect spherical drops which solidified in the 
air and were collected on sheets of paper spread around the sprayer. 

Making glass spheres proved a much harder task. After many trials the following 
procedure was found satisfactory. A thin glass rod was heated to a white heat in a point- 
ed flame and swiftly drawn out so that a fine fibre about 0.02 mm thick was formed. 
The one end of the fibre was introduced into the pointed flame, it melted instantly 
and formed a small sphere hanging by the fibre. It was brought under the microscope 
with the micrometer eyepiece to ascertain its diameter. If it was smaller than required 
the operation was repeated. The sphere was then cut off from the fibre and touched 
with the pointed end of a thin graphite rod that was previously heated in the Bunsen 
flame and kept some time in water. The sphere stuck to graphite and was brought 
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into the pointed flame directed downwards, heated to a white heat, and shaken off 
into a glass beaker. It was then brought again under the microscope for final examina- 
tion. It is worth mentioning that the spheres after striking the bottom of the beaker 
made of a different sort of glass became so strongly electrified that it was difficult to 
remove them from the beaker. 


oo. < 
6 600g @ 80000 


Fig. 6. Wax spheres in the Kundt Fig. 7. Glass spheres in the Kundt 
tube at frequency > = 250 c. p. s. tube at frequency > = 250 c. p. s. 


The spheres a few tenths of a mm in diameter were introduced into the, Kundt 
tube and the air set in vibration. The spheres formed up into lines normal to the 
axis of the tube as shown in Figs. 6 and 7, relating to spheres of 0.5 mm diameter 
made of wax and glass respectively. The frequency was in this case v = 250 c. p. s. 
and the potential difference on the voice 


coil of the loudspeaker U — 65 V. At ч жы a 
this intensity the spheres are distinctly e ооо Oo ы 
separated. 

In order to produce one or two 


rows of spheres at a given place it is © оооо© © 
sufficient to put into the tube two bent | 
strips of lead 3—4 mm thick. They  . а тЫ 
prevent the spheres from wandering 


along the tube and exert a stabilizing о o ооо © ex 


action on the particles within the en- 
closure. 


In some of such experiments the 


Carnauba wax spheres dispersed when ива © P oO о 
the vibrations of the air ceased. Pushing | о ын г оо Тт; © | 
one sphere towards another caused the Anin 5 n 


i 


71 


| 


latter to fly away. Nearing of an electri- Fig. 8. Glass spheres at increasing sound 
fied ebonite rod repelled the spheres intensity. 


which returned to their former places 
when the rod was removed. Obviously electrostatic forces may sometimes appear 


and perhaps influence the formation of ridges. Andrade experimenting with wax 
spheres noticed some irregularities in the ridges. He repeated the experiments with 


Е 
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the tube silvered inside and earthed and came to the conclusion that the influence 
of electrostatic forces was rather doubtful. The question deserves perhaps a closer 
study. 

In order to avoid the disturbing effects of electrostatic forces in the experiments 
with spheres resting on glass the spheres were prepared from the same glass as the 
tube and the air in the tube was ionized by means of a few radium needles used in radio- 
therapy. Fig. 8 shows the row of seven spheres formed up at frequency » — 240 c. p. s. 
and four voltages U = 30, 35, 40 and 45 V on the voice coil of the loudspeaker. 
The increase of the mutual 
repulsion of the spheres with 
the intensity of the sound wave 
is here beyond question. 

At larger distances the 
spheres attract each other. The 
force of attraction, however, is 


= mme EEE Үр 


УМУМЙ 


smaller and is more pronounc- 
ed for particles of a substance 
of low specific gravity, e. g. 
cork or pith. 

To determine approxima- 
tely the repulsive forces acting 
between the spheres in trans- 
verse position the following 


Fig. 9. Experimental arrangement experiments with spheres on . 
for measuring the force of repulsion. the inclined plane were carried 
out. 


A glass plate was attached to the brass ring R fitting closely into the Kundt tube. 
The ring could be rotated about the horizontal axis and the plate inclined at different 
angles to the level. A long pointer was attached to the protruding end of the ring. 
The displacement s of the pointer could be read off on the horizontal mm scale (Fig. 9.). 
The distance from the end of the pointer being 185 mm the angle of inclination was 
found from the equation s = 185 tga. In the centre of the glass plate a small sphere S, 
of 0.33(4) mm diameter was pasted up and another one S, of 0.34(6) mm diameter was 
placed close to it. The spheres, the plate and the rectangular enclosure preventing 
the free sphere from falling from the plate were made of the same glass sort. The 
frequency of the sound wave was у = 253 c. p. s. and the voltage on the voice coil 
of the loudspeaker U — 60 V. When the vibrations obtained in the tube the ring R 
was rotated by small steps and snapshots of the spheres were taken with the photo- 
graphic camera ,Practica^ adapted for large magnifications. The diameters of the 
spheres and their distances were measured on the film by means of a Zeiss comparator. 
True distances were then calculated from the known magnification as the diametres 
of the spheres could be measured with the aid of the same intrument. 
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The results of these measurements are given in Table I. 


TABLEI 


Force of repulsion of spheres resting on glass 


Displacement | Inclination Distance True distance Distance Force 
No | of the pointer | of the plate | of the spheres | of the spheres | of the centres of repulsion 
on the film 

| 5 а r F 

1 0 mm 0.09 3.9 mm 0.936 1.276 

2 2 0.77 1.82 0.437 0.777 0.6х 10 3dyn 

3 5 1.54 1.65 0.396 0.736 1.3 

4 10 3 122 0.293 0.633 2.6 

5 15 4.6 0.76 0.182 0.522 3.8 

6 20 6.2 0.59 0.142 0.482 5a 

7 25 ТТ 0.51 0.122 0.462 6.4 

8 30 92 0.53 0.127 0.467 7.6 

9 35 10.7 0.43 0.103 0.443 8.9 

10 40 12.2 0.35 0.084. 0.424. 10.1 

11 45 13.7 0.25 0.060 0.400 11.3 

12 50 IS 0.14 0.034 0.374 12.4 

18 55 16.5 0 0 0.340 


In the last column the force of repulsion is given. It was calculated from the 

equation 
E= = ғ dg sin а 

expressing the equality of the force of repulsion and the component of the weight 
of the free sphere parallel to the plate, the effect of friction being neglected. The 
density of glass d = 2.5 and the acceleration of gravity g = 981. | 

Table II contains the results of experiments with somewhat larger spheres of 
equal diameters 2r, = 2r, = 0.39(2) mm and higher frequency v = 380 c. p. s., 
the voltage energizing the loudspeaker being as before U — 60 V. The angles of 
inclination a of the plate were in this case determined by the mirror and scale method. 
The mirror was attached to the protruding end of the plate; the scale was fixed in 
vertical position at 1392 mm distance from the mirror. The deflection of the light 
spot when the glass plate was inclined is now s — 1392 tg 2a, a being the angle 
of inclination. | 

As seen from the tables and corresponding graphs (Figs. 10 and 11) the терш 
increases quickly with decreasing distance of the spheres. According to König’s 
theory the spheres should attract each other in transverse position and the force 
between them should vary inversely as the fourth power of the distance between their 
centre. This theory does not apply to vibrations of large amplitude with accompanying 
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TABLE II 
Force of repulsion of spheres resting on glass 
Displacement | Inclination Distance True distance| Distance Force 
No |ofthe pointer | of the plate | of the spheres | of the spheres | of the centres of repulsion 
on the film 
s a r H 
| ы. 26 Е 

1 28 mm 0.58 0.93mm 0.246 mm 0.638 mm | 0.78x 107%4уп 

2 30 0.62 1.22 0.323 0.715 0.83 

3 56 1.15 1.47 0.390 0.782 1.6 

4 58 12 0.86 0.228 0.620 1.6 

5 88 1.8 0.71 0.188 0.580 2.4 

6 90 1.9 0.73 0.193 0.585 2:5 

H 117 2.4 0.62 0.165 0.557 3.2 

8 118 2.4 0.57 0.151 0.543 3.1 

9 146 3.0 0.58 0.154 0,546 4.0 
10 148 3.0 0.43 0.114 0.506 4.1 
11 178 3.6 0.44 0.117 0.509 4.9 
12 179 3.7 0.42 0.111 0.503 4.9 
13 208 42 0.38 0.101 0.493 5.7 
14 209 4.3 0.45 0.119 0.511 5.7 
15 237 4.8 0,34 0.090 0.482 6.5 
16 237 4.8 0.45 0.119 0,511 6.5 
17 268 5.5 0.23 0.061 0.453 7.3 
18 268 5.5 0.28 0.074 0.466 13 
19 298 6.0 0.30 0.079 0.471 8.1 
20 298 6.0 0.30 0.079 0.471 8.1 
21 304 6.2 0.28 0.074 0.466 8.3 
22 328 6.6 0.37 0.098 0.490 8.9 
23 328 6.6 0.27 0.069 0.461 8.9 
24 358 7.2 0.25 0.066 0.458 9.7 
25 358 2 0.25 0.066 0.458 9.7 
26 388 7.8 0.21 0.056 0.448 10.6 
2 388 7.8 0.23 0.061 0.453 10.6 
28 412 8.3 0.06 0.016 0.408 11.1 
29 418 8.4 0.16 0.042 0.434 11.2 
30 438 8.7 0.12 0.032 0.424 11.7 
3l 448 8.9 . 0.18 0.048 0.440 12.0 
32 453 9.0 0.20 0.053 0.445 121 
33 456 9.1 0.07 0.019 0.411 12.2 
34 413 9.4 0.07 0.019 0.411 12.6 
35 413° 9.4 0.00 0.000 0.392 12.6 
36 498 9.8 0.00 0.000 0.392 13.2 


vortices which, as shown by Andrade, are essential for the formation of ridges. 
Nevertheless it seems interesting to try the formula 
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0,010 


Inclination of the plate 
Repulsion in dynes 


| 1 12 mm 
distance 


Fig. 10. The force of repulsion as a function of distance. 
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Pc A 
for the force of repulsion in the above experiments. 
Taking the logarithms we get the equation 
log F = — n log r + loge 


from which n can be found graphically. À | 
From Fig. 12 it follows that plotting the logarithms of force against the loga- 
rithms of distance between the centres of the spheres gives approximately a straight 


logr log r 


-14 -12 -10 


Fig. 12. The logarithms of the force as a function Fig. 13. The logarithms of the force as a function 
of the logarithms of the distance between the cen- of the logarithms of the distance between the | 
tres according to Table I. centres according to Table II. 


line corresponding to n = 4. This applies to experiments listed in Table I. From 
the curve of Fig. 13 relating to experiments with larger spheres (Table IT) follows 
a somewhat larger value of m. 

The repulsion taking place between spheres of different size it seemed probable 
that at sufficient sound intensity a sphere should be repelled by a glass plate. To prove 
this a glass plate was fixed to the bottom of the tube vertically and parallelly to the 
axis. А glass sphere having 0.35 mm in diameter was placed close to the plate so that 
it was seen to touch its virtual image (Fig. 14, a and e). When sufficiently strong 
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vibrations were excited in the tube the sphere moved away from the plate as shown 
in Fig. 14, b, c, d relating to the voltage of 40, 60 and 80 V on the voice coil. 
The above experiments suggest that the walls of the tube may also repel small 
particles in strong acoustic fields, the bottom tending to lift them up. It was found 
impossible to overcome thus the force of gravity acting on glass spheres a few tenths 
of a mm in diameter, which seems plausible in view of the 
smallness of the repelling force exhibited in Tables I and II. 


It should be possible, however, to show the lifting effect of = 
the bottom of the tube on small spheres made of very light d 
material such as cork or sunflower pith. The preparation of 

such spheres being very difficult, experiments were perform- € 
ed with thin disks of 0.7 mm diameter cut from pith. Such ; 

a disk lying on the bottom of the tube stands up transversely b 
to the axis of the tube when a sufficiently strong acoustic " 


field sets in, and is lifted perceptibly in the air when the 
intensity of sound is further increased (Fig. 15). If several ^ Fig, 14. Repulsion by 
disks are introduced into the tube the effect is similar (Fig. a glass plate. 
16). The disks form a row normal to the axis of the tube 

and are plainly separated from one another. It frequently happens that one sphere 
is lifted by two neighbouring ones the higher the larger the amplitude of the air 
vibrations. 


Fig. 15. Lifting of a disk Fig. 16. Lifting of three disks in the Kundt tube. 
in. the Kundt tube. 


Similar effects were observed on small spheres of the substance used by dentists 
for making artificial jaws. 

It should be noticed that neither single particles nor groups of particles stand 
still in the tube when sound waves obtain, but wander to and fro along the tube. They 
can be kept quiet by nearing from above a glass rod bent downwards. 

III Forces between spheres suspended in the Kundt tube. 

In these experiménts one sphere $; was suspended on a fine glass fibre attached 
to a thin brass ring fitting tightly into the tube. Another sphere S, was fixed on a thick 


TM 
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glass fibre connected with the micromanipulator M and could be moved with great 
precision in any direction. When the air in the tube is set in vibrations the thin glass 
fibre starts vibrating with the same frequency displacing somewhat the sphere S, 
from its position of equilibrium to a new one which may.be conveniently called the 
„dynamical zero“. If the spheres are of the same or similar size as those employed 2 
the experiments referred to above (Tables I and П), and the sound intensity is suffi- 
ciently large, the repulsion of the suspended sphere 5, under the influence of the 
sphere S, approaching it transversally may be easily observed either visually or photo- 
graphically. The dynamical zero being however not sufficiently stable, exact measu- 
rements cannot be made in this way. To stabilize the dynamical zero the following 
device was used. 


mm 


02 


> 
displacement 


repulsion 


distance 


III III III А N 
а 


d 
attraction 


Fig. 17. The ring supporting the Fig. 18. Displacement of the suspended sphere as a 
glass fibre. function of distance between spheres having diameters 
0.34 mm and 0.35 mm. 


To the ring supporting the glass thread (Fig. 17) was sealed a short copper wire 
C 2 mm in diameter on which a glass tube G 20 mm long, 3.5 mm internal diameter 
was fixed by means of a thin rubber collar. The fine glass fibre supporting the sphere 5; 
was stuck to the copper wire with picein and the tube was filled with glycerine to the 
brim. It could be placed in the tube upright or upside down, the surface tension of 
glycerine preventing its being spilled. In the former position the sensitiveness of 
the suspension was greatér the gravity tending to counterbalance the restoring force 
of the glass fibre. 

In spite of these precautions the dynamical zero was not very stable for spheres 
‚ about one third mm in diameter. This entailed an error in the determination of the 
displacement of the movable sphere S, caused by the presence of the fixed sphere S, 


=e E 4 
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In order to eliminate as far as possible the disturbing effect of the unstable zero nine 


photograhps were employed to find the distance between the spheres and the dis- 
placement of the suspended sphere, viz. 


1, 2, 3 in the absence of S, 
4, 5, 6 „ , presence „ 
7,8,9 „ , absence 


The results of measurements with glass spheres having the diameters 2n = 
— 0.34 mm and 2r, — 0.35 mm in the sound waves of frequency — 375 c. p. 5. 
and at 60 V potential difference on the voice coil of the loudspeaker are shown 
in Fig. 18, where the displacement of the suspended sphere is plotted against the 
distance between the spheres. 

The displacement of the sphere 5; from the dynamical zero may be taken for the 
approximate measure of the force of repulsion exerted on it by S. Displacements 
directed away from S, are considered positive (repulsion), those directed towards 
it — negative (attraction). 


29 


The measurements display considerable scattering, nevertheless the graph shows 
unmistakably the repulsion of spheres at small distances and the attraction at larger 
ones. The attraction just perceptible at distances exceeding five times the diameter of 
the spheres increases with decreasing distance, reaches a maximum and vanishes at 
a distance of about 0.3 of the diameter to turn into repulsion growing rapidly with 
further approach of the spheres. 4 

The above experiments show that under the specified conditions Cook’s assertion 
concerning the force acting on spheres placed transversely to the vibration vector 
is undoubtedly right. On the other hand Andrade’s investigations can hardly be put 
in doubt. The experiments described below explain the reason of this discrepancy. 
They show that the force in question depends not only on the relative position of 
the spheres but on their size as well, and on the frequency and intensity of the sound 
waves. Now, Andrade investigated the behaviour of suspended cylinders?) of 30 mm 
in length and 1.5 mm in diameter at a sound frequency of 512 c. p. s., whereas Cook 
experimented with spheres having diameters ranging from 0.2 mm to 0.9 mm at V = 
= 1287. The relative intensity of sound is not known. No wonder that the results of 
experiments under so widely different conditions were contradictory. 

The influence of the size of the spheres, as well as of the intensity and frequency 
of the air vibrations was investigated employing somewhat larger spheres than in 
the previous experiments, suspending them on thicker glass fibres and stabilizing the 
dynamical zero by means of more viscous fluids (mixtures of honey and glycerine) 
surrounding the fibre. Greater accuracy was thus obtained in shorter time for the 
positions of the spheres could be determined visually by means of the long focus 
microscope with filar micrometer eye-piece. 
mun pem SIT ee Te ee 

?) The experiments with spheres are only just mentioned by Andrade. 
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1) Different spheres in the same sound field. 
The behaviour of spheres having diameters from 1.04 mm to 5.45 mm was 


observed in the same sound field at V = 384 c. p. s. and U = 60 V. The smaller 
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Fig. 19. Displacement of the suspended sphere 
(2r,= 1.03 mm) as a function of distance from 
the fixed sphere (2r, = 1.32 mm). 


distance 


mm 
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Fig. 21. Displacement of the suspended 

sphere (2r, = 1.55 mm) as a function of 

distance from the fixed sphere (2r, — 1.32 
mm). 
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Fgi. 20. Displacement of the suspended sphere 
(2r,— 1.17 mm) as a function of distance from 
the fixed sphere (27, = 1,32 mm). 
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е 
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Fig. 22. Displacement of the suspended 
sphere (2r, — 1.92 mm) caused by fixed 
spheres of different size (2r, — 1.32, 2.31, 
3.5, 5.45 mm) as a function of distance. 


suspended spheres were made of picein, the larger ones of alder wood. The fixed sphere 
was made of glass. The graphs in Figs. 19, 20, 21 show the results of measure- 
ments with three pairs of spheres having the diameters 1.04 mm and 1.32 mm, 1.17 mm 
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and 1.32 mm, 1.55 and 1.32 mm. The scale of distances is the same in the above graphs, 
but that of the displacements is different. To equal displacements in different graphs 
do not correspond equal forces because the thickness of the glass fibre and the weight 
of the suspended sphere are different. The forces acling between different pairs of 
spheres are therefore not directly comparable. 

All curves display one characteristic minimum representing the maximum of 
attraction at a distance smaller than half a diameter of the spheres. The curves rise 
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Fig. 23. Influence of the intensity of air Fig. 24. Infiuence of the frequency of the 
vibrations on the force acting between two sound field on the force acting between the 
spheres. spheres. 


slowly for larger distances and faster for smaller ones. In Figs. 19 and 20 relating to 
smaller suspended spheres the curve crosses the X-axis, which means that at small 
distances the force becomes one of repulsion, just as in the case illustrated in Fig. 18. 
In Fig. 21 however corresponding to a larger suspended sphere (2r, = 1.55 mm) 
the curve has as it were no chance of crossing the X-axis before the spheres touch. 
This is an analogous case to that of Andrade's ivory cylinders. · 

In Fig. 22 are given the displacements of the wooden sphere of 1.92 mm diameter 
caused by four glass spheres having the diameters 1.32, 2.31, 3.5, and 5.45 mm. 
The wooden sphere was suspended on the same glass fibre in all experiments, so 
that its displacements may serve as a measure of force in arbitrary units. It is seen 
from the graphs that the attraction increases at all distances with increasing size of 
the spheres. It does so in particular at the contact of the spheres. In view ofthe con- 
vention of signs made above we may say: The force at the contact of two spheres 
situated transversely to the vibration vector decreases with increasing size of the 


spheres. 
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2) Influence of the intensity of sound. 

Two spheres having diameters of 1.56 mm and 1.50 mm were observed at constant 
frequency v = 380 c.p.s. and two intensities of sound corresponding to 60 and 90 У on the 
voice coil of the loudspeaker. The results given in Fig. 23 show that increasing of 
sound intensity entails a diminution in the attractive force which may even change 
into a repulsive one at small distances between the spheres. Hence, the force at the 
contact of two spheres situated transversely to the vibration vector increases with 
the intensity of air vibrations. 

3) Influence of the frequency of sound. 

The behaviour of the above two spheres at vibration frequencies 184, 380, 553 
с. p. s. and at constant voltage U = 60 V on the voice coil is shown in Fig. 24. 
It follows from these graphs that the effect of increasing frequency is similar to that 
of increasing the size ofthe spheres. Hence increasing of frequency of air vibrations, 
the voltage on the voice coil being constant, diminishes the force at the contact of 
two spheres situated transversely to the vibration vector. 

Similar results were obtained if the current energizing the loudspeaker was 
adjusted so as to keep constant the voltage induced in the microphone placed close 
to the open end of the Kundt tube. Hence taking into account a fairly constant sensi- 
tivity of the microphone used it may be said that increasing frequency at constant 
intensity of the sound field diminishes the force between the spheres with their line of 
centres transverse to the vibration vector. 

The author wishes to thank Professor S. Ziemecki for his valuable aid in making 
this work possible. Special thanks are due to Mr. B. Adamczyk for his untiring assistance 
in planning and carrying out the experiments. 


KPATKOE СОДЕРЖАНИЕ 


B. Сташевский, О взаимодействии двух шариков B вибрирующем воздухе. 


Было исследовано взаимодействие двух шариков помещенных в трубке 
Кундта поперек трубки. Обнаружено, что сила, действующия между шарика- 
ми, может быть равно притягательная, как и отталкивающая, в зависимости OT 
диаметров шариков, от частоты и от напряженности звукового поля. Этим объяс- 
няется противоречие результатов опытов Кука и Андраде, осуществленных в да- 
леко не одинаковых условиях. 
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LETTERS TO THE EDITOR 


THE DEUTERON POLARIZABILITY AND ITS EFFECT 
ON THE RUTHERFORD SCATTERING 


By JERZY SAWICKI 


Institute of Theoretical Physics, Warsaw 


(received January 25 1954) 


Kruse, Malenka and Ramsey (1953) have published two papers in which the 
deuteron polarizability in an electrical field and the deviations from the Rutherford 
scattering due to the deuteron spatial structure were discussed, and the cross-section 
was evaluated by classical method, the deuteron energy in the Coulomb field of the 
nucleus being ; 


aZ?e? 


Де? 
m Ufa eem E is . 10-39cm3 
Е = Mv? + R + у (К), where 7 (К) JRI (a = 0,56 ° 107*9?cm?) 
is the energy correction due to polarization. The classical theory of scattering given 
Ze? ZC 1 
by these authors is restricted to the range: л = STE aE 1. The aim of the 


present work is to obtain 7 (R) with the use of another method and to discuss 


quantum-mechanically the scattering deviations. 
Let us start from the French and Goldberger (1952) equation: 


h Le E E a 1 HE 
КЕ DE, =: R т МА vo | Ф (К, 7) = Ze? = EI Ф (К, т), 
2 (1) 


10 t n Team 7 — radius vectors of the proton and the neutron), 
E, is the deuteron full energy, U (г) — the n—p interaction potential and M the 
nucleon mass. On expanding Ф (R, г) in a set of роо eigenfunctions gg (В, 7) 


for the internal n—p motion: 


D (R, = V, (R) y (Ror) + f dE Vx (R) v (Re? @ 
(225) 


where Р = 
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we get for the ground state wave function 


2 е. СЕ 
Ес + U(r) —e + (к) |») = 
A 2 
Ze? Гран x] Фо(К,т) = Е (К) (3) 
mo 


If x, (7) is the deuteron ground-state wave function, then by Kirkwood's variational 


method (see Hellmann 1937) we get 


1 (R) = Уы ree P ere" 
@ 


7710 lo m P (№ 
5a  Ze?rcosO á 
where: V (К, г) әк ОИ) S45 [D(V loo = < 0|(у„Ё) |0>. 
wa | | 
Here Г» = 0 and — 7 (R) = ——„————°9——. To compare the results obtained by 


h2 
17 Dlo 
our method and by the K. M. R. method, we use the Hulthen function employed by 
these authors: 


© =~ (ere) 


Where û =O LZ Қа, I TIRES = 0,675: OSE Cm 
aZ?e? 

Thus we get — 7 (R) = 3g ^ 
2 eM B?r? (B -- D 1 4 ICE 
uj E38 = - 10739 3 
Dah? ДА a (t D А tee 
as compared with the K. М. В. result а = 0,56 - 10-39 стз, For the Wilson function 

D = Jor BET ze - 10-13 з еМ ба 
Xo(7) ОЕ; with o^! = 4,5 · 10-13 cm, we get a A? 79 

The exact computations based on the method of Kirkwood being equivalent 
to the second order perturbation theory seem to give a better result for a than the 
K. M. R. second order perturbation calculations with replacing of the deuteron wave 
equation Green's function by the free particle Green's function. On inserting „adia- 
batically^ the term with 7 (R) into the equation for Y^ (В), we get for У (Б) 


where а = 


= 0,21 - 10-39 cm. 


h? Ze? = os = 
[erm е |6 =n] > ru 6 
Eq. (5) is valid for R >a, а = 1,5 41'3. 10-1? cm being the radius of the nucleus 
(the Coulomb barrier) and y, (R) the Coulomb wave function. 
On employing the asymptotical form of Green’s function for large R, we have 
for the solution of the inhomogeneous equation 
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~ 


1 5 i(kR — nIn2kR) 
y pico Um AU — = ж /— 
0 da R 3% 7 (r) v1 ©) va (F) dr, (6) 


where 7] (г) = 2 (г) for r Еа) 0 form ча („eutting о 
the barrier), y, (7) = N eser, Е, (— in, 1; 51 (ki r— kı 7)), 2 (r) = N e-ikir 17 (— in, 
1, i (kar + ka 7)), N? =e" I (1 + in)?, kı = F, ka = 23 k = 28 ‚МБ = 
= Kk? cos О, © being the scattering angle, „F; — the confluent hypergeometrical func- 
tion. The evaluation of the last integral is extremely difficult. However, in order to 


gain some information concerning the lower limits of the deviation 6 = ie x] 


б 
к 
for some cases, the character of the dependence of 6 оп n and the incident deuteron 


energy, the angular distribution — we can evaluate (6) on replacing 7] (т) by у (r + a). 

1 
(r + а)“ ln 
sequence of integration. We perform the 7 Are. on applying the Sommer- 
feld method generalized for the case 4 + 0. We finally obtain the scattering amplitude 


We apply the Laplace transform 


q? e “e "dq and interchange the 


AO == E. aZ?ete-"^ Г (1 + in)? = I (Ө) where 
a 2 2. 19k 2in 
I(0) - [== (4-1) 7 ч 
б G + 4k? sin? 2) 


| Ae 
F(- in; —in, 1; а sin? | adq 
or with the aid оЁ а well known formula 


F(—ngy = па, l; —x) = (1-р ола Е(1--1-- пә, 1; —2): 


2in 2‏ م 
I(0) = fe С — ) (1-4 zm (1 + in, 1+ in, l; un) adq‏ 
б‏ 


F being the hypergeometrical function. If f, (O) is the „pure“ Coulumb scattering 
amplitude, we conclude, that the ratio of the differential cross-section do to the purely 
Rutherford cross-section dog is 


(7) 


йо = 2. = ГӨ) "s -лп E а 2 k2 2 EE ton 21 0 

dos = 1 а AE A FLO = Za E e е sin ( ) (9) 
From the above we see that A = 0 for © = 0. For small scattering B A is negligible. 

The approximate form of A was obtained for large scattering angles and Z>1 
with the use of two different methods: one for n< 2,5 and the other for considerably 
le ndun Жыслал АЙ з ANA Bee шз. a „ш Sof fae 


1 The main results are the same for some other forms of the polarization potential. 
ri 
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larger n. In both cases after approximate integrations we obtained A as weakly depen- 


em 


dent on л (the factor does not appear in A) .E. g. for n< 2,5 (large energies) 


shan 
with use of the analytical continuation of F (1 + in, 1 + in, 1; 2 for | z | > 1 we 
find for © = 180° with a good approximation: 


zonae 2 ат 
әгі E 2 (nay — 49 (3) + 3y (2) + In 7 —ce— y (in) «zl хес 


у = kn, y (z) is the „ря“ Euler function and С - Euler constant. For а = 
= 1,5 ?/938 - 10-13 cm, Z = 92, n = 2, а = 0,32 - 10-9? стз, we gotö = — 0,174% 


for Ө = 180? (maximum of | 0|), д = — 0,172% for Ө = 140°; û is negative 
because the attractive potential 7 (R) causes a decrease of backward scattering. 


Ake 


For considerably larger n (smaller energies) 6 was evaluated with use of the form 
of F (—in, —in, 1, z,) obtained after Sommerfeld with use of the ,,steepest descent“ 
method. For the above values of а, Z, a and п = 6,3 (energy ~ 10 MeV, the К. М. В. 
case), Ө = 180°, we get: д = —0,09% and ô = —0,15% for а = 0,56 ° 10°3’cm? 
(the K. M. R. value). The last values of 6 are only slightly different from those for 
n = 2. Our last values’of (—6) are much smaller than that obtained from the К. М. К. 
classical theory (—6 = 4,2%). In reality the accurate values of (—6) are considerably 
larger than those obtained here because of taking the potential V, = — 7 (r + а) 
instead of V, = — 7 (г) in (6)!). We observe, that f „йт, the „volume“ of V, was 
thus decreased and thus smaller values of the Coulomb wave functions were used, 
with r instead of r -+ a in their arguments in the most important area of integration 
(E d): 

For very small energies an approximate „equivalence“ of 7 (г) and 7 (r + a’), 
a’ being the „equivalent radius“, may be obtained with a’ = a/3 (equal „volumes“ 
and corresponding values of д for 7 (г) and y (r -- a). If we applied it for the larger 
energy in our last case (ô = --0,15%), we would find ô = —1,3%). 

I wish to express my sincere thanks to Dr R. Kolodziejski for suggesting the 
problem, and many valuable discussions. I am also very indebted to Professor 

‚1. Infeld and Dr J. Werle for helpful criticism. 
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